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Abstract 

Galactic globular clusters are old, dense star systems typically containing 10*-10^ stars. 
As an old population of stars, globular clusters contain many collapsed and degenerate ob- 
jects. As a dense population of stars, globular clusters are the scene of many interesting close 
dynamical interactions between stars. These dynamical interactions can alter the evolution of 
individual stars and can produce tight binary systems containing one or two compact objects. 
In this review, we discuss theoretical models of globular cluster evolution and binary evolu- 
tion, techniques for simulating this evolution that leads to relativistic binaries, and current and 
possible future observational evidence for this population. Our discussion of globular cluster 
evolution will focus on the processes that boost the production of hard binary systems and 
the subsequent interaction of these binaries that can alter the properties of both bodies and 
can lead to exotic objects. Direct A''-body integrations and Fokker-Planck simulations of the 
evolution of globular clusters that incorporate tidal interactions and lead to predictions of rel- 
ativistic binary populations are also discussed. We discuss the current observational evidence 
for cataclysmic variables, millisecond pulsars, and low-mass X-ray binaries as well as possible 
future detection of relativistic binaries with gravitational radiation. 
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1 Introduction 



For the purposes of this review, relativistic binaries are systems containing two stellar mass de- 
generate or collapsed objects that are in close orbits. In the Galactic field, where these systems 
evolve in relative isolation, their final properties are set solely by their initial conditions and are 
the result of mass transfer, common envelope evolution, or other interactions that may interrupt 
the course of single star evolution due to the presence of a nearby neighbor. When considered as a 
fraction of the total stellar mass, the number of relativistic binaries in Galactic globular clusters is 
overrepresented compared to the Galactic field. Thus, the dynamical interactions found in the en- 
vironment of dense stellar clusters provide additional channels for the formation of these systems. 
Relativistic binaries reveal themselves observationally as UV or X-ray sources and are potential 
sources of gravitational radiation. 

This review will concentrate on the Galactic globular cluster system for the bulk of the text. 
We shall touch on extra-galactic globular cluster systems briefiy in the sections on observations 
(Section [3]) and gravitational radiation prospects (Section 15]). 

We begin in Section [2] by looking at the physical structure and general history of the galactic 
globular cluster system that leads to the concentration of evolved stars, stellar remnants, and binary 
systems in the cores of these clusters. Current observations of globular clusters that have revealed 
numerous populations of relativistic binaries and their tracers are presented in Section [3l We 
also consider the prospects for future observations in this rapidly changing area. Many relativistic 
binaries are the product of stellar evolution in close binaries. In Section|4l we will look at how mass 
transfer between one star and a nearby companion can dramatically alter the evolution of both 
stars. The enhanced production of relativistic binaries in globular clusters results from dynamical 
processes that drive binaries toward tighter orbits and that preferentially exchange more massive 
and degenerate objects into binary systems. Numerical simulations of globular cluster evolution, 
which can be used to predict the rate at which relativistic binaries are formed, are discussed 
in Section [5j These models can be compared with the observable members of the population of 
relativistic binaries in order to try and constrain the entire population. Finally, we conclude with a 
brief discussion of the prospects for observing these systems in gravitational radiation in Section El 

Readers interested in further studies of the structure and evolution of globular clusters are in- 
vited to look at Binney and Tremaine |40j , Spitzer [356] , and Volumes I and II of Padmanabhan's 
Theoretical Astrophysics [28511286] for a good introduction to the physical processes involved. Re- 
view articles of Meylan and Heggie |257) and Meylan [256] also provide a comprehensive look at 
the internal dynamics of globular clusters. Although our focus is mainly on the Galactic globular 
cluster system, the physics of globular cluster systems associated with other galaxies is well covered 
in the review article by Harris [1501 as well as his lecture notes from the Saas-Fee course on star 
clusters [48]. Carney has a thorough introduction to evolution of stars in globular clusters [49] . 
An observational perspective on the role of binaries in globular clusters is presented in an excel- 
lent review by Bailyn [16] . while a good introduction to the details of observing binary systems 
in general can be found in An Introduction to Close Binary Stars [172] . Although slightly out of 
date, the review of binaries in globular clusters by Hut et al. [194] is an excellent introduction to 
the interaction between globular cluster dynamics and binary evolution, as is a short article on 
globular cluster binaries by McMillan, Pryor, and Phinney [254] . Rappaport et al. [317] and Rasio 
et al. [319] have written reviews of numerical simulations of binary populations in globular clus- 
ters. An excellent introduction to the astrophysics and numerical techniques relevant to globular 
cluster dynamics can be found in the book by Heggie and Hut [155] . Finally, a shorter and more 
observationally focused review of compact objects in globular clusters can be found in Maccarone 
and Knigge [244] . 
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2 Globular Clusters 



Globular clusters are gravitationally bound associations of 10'' — 10'' stars, distinct both from their 
smaller cousins, open clusters, and the larger, dark matter dominated dwarf galaxies that populate 
the low-mass end of the cosmological web of structure. Globular clusters are normally associated 
with a host galaxy and most galaxies, including the Milky Way, are surrounded and penetrated by 
a globular cluster system. A good estimate of the number of globular clusters in the Milky Way 
is the frequently updated catalogue by Harris jl51] . which has 157 entries as of 2010. Although 
fairly complete, a few new clusters have been discovered in recent years at low Galactic latitudes 
[1871 1220) and there may be more hidden behind the galactic disc and bulge. The distribution of 
known globular clusters in the Galaxy is given in Figure [1] Other galaxies contain many more 
globular clusters (M87 alone may have over 10 000 11521 ) and the luminosity of globular clusters 
normalized to the host galaxy luminosity seems to be higher for elliptical than disc galaxies [44) . 




Figure 1: Globular cluster distribution about the galaxy. Positions are from Harris [151] and are 
plotted as black circles on top of the COBE FIRAS 2.2 micron map of the Galaxy using a Mollweide 
projection. Figure taken from Brian Chaboyer's website [52 . 

Milky way globular clusters are old, having typical ages of 13 Gyr and an age spread of less 
than 5 Gyr [50] . This is on the order of the age of the Galaxy itself, thus Galactic globular clusters 
are thought to be left over from its formation. By contrast other galaxies such as the small and 
large Magellanic clouds (SMC and LMC) have intermediate age globular clusters (< 3 Gyr old, 
e.g. |2741 1268] ) and in some galaxy mergers, such as the Antennae, massive star- forming regions 
that may become globular clusters are observed ^3). Taken together, this implies that globular 
clusters of all ages are relatively common objects in the universe. 

2.1 Stellar Populations in Globular Clusters 

Most of the detailed information on stellar populations in globular clusters comes from those in 
the Milky Way since only they are close enough for stars to be individually resolved. The stars 
in individual Galactic globular clusters all tend to have the same iron content |137J so globular 
clusters are thought to be internally chemically homogeneous. The colour-magnitude diagram 
(CMDs) for most Galactic globular clusters (e.g. M80, Figure [2]) also indicate a single stellar 
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population with a distinct main-sequence, main-sequence turn-off, horizontal and giant branch. 
The single main sequence turn-off in particular indicates a co-eval stellar population. This leads 
to a so-called "simple stellar population" model for globular clusters where all stars have the same 
composition and age and differ only by their masses, set by the initial mass function (IMF). This 
simple picture has been challenged in recent years as observations have shown systematic star- 
to-star light element variations in globular clusters [137) . This has lead to the development of 
self-enrichment models where a single globular cluster experiences several bursts of star formation, 
each enriched by pollution from the previous generation . How multiple populations will affect 
the CMD of a globular cluster is shown in Figure[3]The importance of these scenarios for relativistic 
binaries has not yet been explored. Here we outline the consequences of the co-eval model. 
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Figure 2: Colour-magnitude diagram for M80. Figure taken from the catalogue of 52 globular 
clusters |328) . The entire catalogue is available at the Padova Globular Cluster Group website [130j . 



The IMF is thought to be universal ;27, and is usually taken to be a power-law of the form 

dN 



dM 



(1) 



where ai can take different values for different mass ranges. For values above ~ IM© Ui is usually 
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Figure 3: CMD for NGC288 showing evidence of two populations, and models incorporating 
different metallicities, helium fraction, and age. The best fit model has an age difference of A Age — 
1.5 Gyr. Figure taken from Roh et al. [327) . 
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assumed to have a single value, the so-called Salpeter slope, of ^ 2.35 [331] . There is much more 
debate about the value of in low-mass regime. One popular choice, the Kroupa IMF, is a broken 
power law with exponent values of 

ao = 0.3 ± 0.7, 0.01 < M/Mq < 0.08 

ai = 1.3 ±0.5, 0.08 < M/Mq < 0.50 , , 

a2 = 2.3 ±0.3, 0.50 < M/Mq < 1.00 ^ ' 

aa = 2.3 ±0.7, 1.00 < Af/Mg 

[2251 1226] ■ Another possibility, introduced by Chabrier [53], uses a log-normal distribution of 
masses below IMq. In all cases it is clear that the IMF strongly favours low masses so stars 
massive enough to form neutron stars (NSs) and black holes (BHs) will be rare. 

Since Galactic globular clusters are old, their stellar populations are highly evolved. Most of 
the stars above ~ O.SMq have already moved off the main-sequence. Those just above O.SMq will 
be on the red giant branch and are readily visible in most optical images of globular clusters, such 
as the image of M80 shown in Figure ]?] Those that are more massive have evolved past the giant 
and horizontal branches and become low-luminosity compact remnants. The fact that the original 
high-mass population is no longer visible produces an intrinsic uncertainty in our knowledge of the 
high-mass end of the IMF of globular cluster stars. 




Figure 4: Hubble Space Telescope photograph of the dense globular cluster M80 (NGC 6093). 
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The picture of Galactic globular cluster stellar populations that emerges from this analysis is of 
a simple, co-eval, chemically homogeneous set of luminous, low-mass population II stars combined 
with a low-luminosity population of high-mass stellar remnants. It is interactions with members 
of this remnant population that will be of particular interest for producing relativistic binaries. 



2.2 The Structure of Globular Clusters 

Globular clusters are roughly spherical N-body systems that can be described by a core-halo 
structure. The core is highly concentrated, reaching densities of up to lO^Mo/pc'^, and strongly 
self-gravitating. The surrounding halo is of much lower density and is less strongly self-gravitating. 
The structure of a globular cluster can be classified using three basic radii: the core radius (r^), 
the half-mass radius (r/i), and the tidal radius (rj). The core radius is often taken to be the radius 
at which the space density drops to about one third, and thus the surface density to one half, of its 
central value |155] (other definitions are possible) . Observationally this corresponds to the radius 
at which the surface brightness drops to half of its central value. The half-mass radius is simply 
the radius that contains half of the mass of the system. The corresponding observational value is 
the half-light radius, which contains half the light of the system (the two radii do not necessarily 
agree). The tidal radius is the radius at which the gravitational field of the host galaxy becomes 
more important than the self gravity of the star cluster. For a given cluster, cluster orbit, and 
galaxy model the tidal radius of the cluster can, in principle, be clearly defined by comparing the 
effect of the galactic versus globular cluster gravitational field on a test mass. Observationally 
tidal radii can be difficult to determine due to the low stellar density of globular cluster halos and 
an imperfect knowledge of the gravitational potential of the host galaxy. Median values for Tc, r^,, 
and rt in the Galaxy are ^ Ipc, ^ 3pc and ~ 35pc respectively [155]. 

There are also two important timescales that characterize globular cluster evolution: the cross- 
ing time (tcr) and the relaxation time (irix)- The crossing time is simply the time required for a star 
traveling at a typical velocity to cross some characteristic cluster radius. Thus tcr ~ R/v where, 
for example, rc or might be typical radii of interest and v could be the velocity dispersion (~ 10 
km s~^). tcr is also, roughly speaking, the orbital timescale for the cluster. For typical values of 
r/i and v, icr for Galactic globular clusters is on the order of 0.1 — IMyr but is longer at the tidal 
radius and much shorter in the core. 

The relaxation time describes how long it takes for orbits to be significantly altered by stellar 
encounters. In particular t^ix is often defined as the time necessary for the velocity of a star to 
change by an order of itself [40] . This can be thought of as the time necessary for a cluster to lose 
the memory of its initial conditions or, more exactly, the time necessary for stellar encounters to 
transform an arbitrary velocity distribution to a Maxwellian |356j . The relaxation time is related 
to the number and strength of encounters and thus to the number density and energy of a typical 
star in the cluster. It can be shown that the mean relaxation time in a globular cluster is given by 
[4011356] 

N 

tiiK oc - — —ta- (3) 
In iV 

Taking a typical value of iV = 10*^ and tcr as before, typical values of trix are 0.1 — IGyr. Thus, 
with ages typically greater than 10 Gyr, Galactic globular clusters are expected to be dynamically 
relaxed objects. In reality the value of irix varies significantly within a globular cluster due to the 
highly inhomogeneous density distribution of the core-halo structure. It is possible for the core of 
a globular cluster to be fully relaxed while the halo remains unrelaxed after 13 Gyr. By making 
various approximations for 1^ it is possible to relate ^rix to local cluster properties. For instance 
the criterion of Meylan and Heggie |2571 1231] : 

^^^^ ^ 0.065(t;^)^/^ 
p{m) ln{jN) 
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relates trix to the local mass density (p), the mass- weighted mean squared velocity ((u^)), and the 
average mass {{m)). The criterion by Spitzer [356] : 

0.138A^i/V3/2 

" (771)1/201/2 in(77V) 

relates trix to a characteristic radius (r), the average mass, and the number of stars in the system. 
In practise rh is normally used for the characteristic radius in the Spitzer criterion and ^rix is 
renamed the half-mass relaxation time, trh- The factor ln7A^ that appears in both definitions 
is called the Coulomb logarithm and describes the relative effectiveness of small and large angle 
collisions. The exact value of 7 is a matter of some debate and ranges from 0.02 — 0.4 depending 
on the mass distribution of the system |124| 1123) . Both criteria are used extensively in stellar 
dynamics. 

On timescales shorter than a relaxation time individual encounters do not govern the overall 
evolution of a stellar system and the granularity of the gravitational potential can be ignored. On 
these timescales the background structure of the cluster can be modelled using a static distribution 
function, /, that describes the probability of finding a star at a particular location in a 6-D position- 
velocity phase space. Formally / depends on position, velocity, mass, and time, thus we have 
f{x,v,m,t). For times less than tiix, however, the evolution of / is described by the coUisionless 
Boltzmann equation: 

|+.-V/-W-| = (6) 

and the explicit time dependence can be removed: df /dt = 0. The gravitational potential, 0, is 
given by Poisson's equation: 

V^(t> = AnGp (7) 

and can be calculated at any position by integrating the distribution function over mass and 
velocity: 

V^</> = AttG I f{x, V, m)(fvdm. (8) 



Solutions to Equation |6] are often described in terms of the relative energy per unit mass, £ = 
^ — 12 where ^' = — (/) -f (/)o is the relative potential and (/)o is defined such that no star has 
an energy less that zero (/ > for > and / = for f < 0). A simple class of solutions to 
Equation [6] are Plummer models |300J : 

/(£) = FZ'I\ (9) 

the stellar-dynamical equivalent of an 77 = 5 ploytrope [155] . Another class of models that admit 
anisotropy and a distribution in angular momentum, L, are known as King-Mitchie models |215l 
1258) . The basic distribution function is given by: 

/(£, L) = pi(27ra2)-3/2 exp(^-^) [e^/^' ~l],£>0, (10) 

where a is the velocity dispersion, ra the anisotropy radius where the velocity distribution changes 
from nearly isotropic to nearly radial, and pi is a constant related to the density. Although not 
as well theoretically supported as the single-mass case, King-Mitchie models have been extended 
to include a spectrum of stellar masses |147) and even external gravitational field |157) . Multi- 
mass King models in particular are often fit to observed globular cluster cluster surface brightness 
profiles in order to determine their masses. A good example of the construction of a multi-mass 
King-Mitchie model is found in the appendix of Miocchi [269) . 
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2.3 The Dynamical Evolution of Globular Clusters 

Although static models can be used to describe the instantaneous structure of a globular cluster 
there is no stable equilibrium for self-gravitating systems [40] and therefore their structure changes 
dramatically over time. Accessible descriptions of globular cluster evolution are given in Hut et 
al. [194] , Meylan and Heggie |257| and Meylan [256] and have also been the subject of several texts 
(e.g. Spitzer |356) and Heggie and Hut |155| ). Here we merely outline some of the more important 
aspects of globular cluster evolution. 

The initial conditions of globular clusters are not well constrained since they seem to form only 
very early in the history of galaxy formation or in major mergers |1181 144] - both situations quite 
different from the environment of our Galaxy today. However it is possible to make some general 
statements. Like all stars, the stars in globular clusters collapse out of molecular gas. Due to their 
small age spread and chemical homogeneity, it seems certain that all the stars in a single globular 
cluster formed out of the same material over a short time period. Thus it seems likely that each 
globular cluster formed by the collapse of a single giant molecular cloud I44j . The details of this 
collapse are unclear. Sub-clumps in the cloud may have collapsed individually and then merged 
later, affecting the course of local star formation. The star formation may have been sufficiently 
extended that early-forming stars could pollute late-forming stars, explaining some observed light- 
element abundance anomalies in otherwise homogeneous star clusters |137j . In particular, the 
star formation efficiency (the fraction of gas converted into stars) is a major issue for globular 
cluster formation. If anything less than 100% of the primordial gas forms stars then the resulting 
globular cluster will be less massive than its parent cloud. This remaining gas will be expelled 
by a combination of radiation pressure from young stars and energy injection from supernovae. 
If the star formation efficiency is low, then the amount of mass loss through gas expulsion can 
leave the cluster out of virial equilibrium and may lead to its immediate dissolution |135l 11361 1384) , 
a process called "infant mortality". It is estimated that over 50% of young clusters in the local 
universe are destroyed in this manner |136j and most surviving clusters will lose a large fraction of 
their stars. Even if the cluster survives the gas expulsion, the rapid change in potential will cause 
the energy of individual stars to change in a mass-independent manner '40'. This process, called 
violent relaxation, means that the positions, velocities, and masses of cluster stars will be initially 
uncorrelated. 

Equipartition of energy dictates that the most massive stars should have the lowest kinetic 
energies [40 ll356j . Thus, as soon as the residual gas is expelled from a young globular cluster, mas- 
sive stars with large kinetic energies will start to transfer this energy to low-mass stars through 
stellar encounters. As the massive stars lose kinetic energy they will sink to the centre of the 
cluster while the low-mass stars gain kinetic energy and move to the halo. This process is known 
as mass-segregation and proceeds on a timescale tms oc mi/{m) j355| 1380"] 1209) . Due to mass- loss 
from stellar evolution, compact remnants rapidly become the most massive objects in globular 
clusters as the star population ages. Thus, they are strongly affected by mass segregation and are 
particularly likely to be found in cluster cores. Mass segregation continues until energy equiparti- 
tion has been achieved. There are, however, initial conditions for which it is impossible to achieve 
energy equipartition j380j 12091 and it is formally impossible to halt mass segregation, leading to 
a singularity. This phenomenon was first noted by Spitzer [355 and is thus called the "Spitzer 
instability" . In reality, the massive objects in such systems form a strongly interacting subsystem, 
dynamically decoupled from the rest of the cluster. Due to their high masses, black holes in star 
clusters are particularly likely to experience the Spitzer instability and this has been the starting 
point for several investigations of BH binaries in star clusters (e.g. |283) ). 

The longer-term evolution of star clusters is driven by two-body relaxation, where the orbits 
of stars are perturbed by encounters with their neighbours. The theory of two-body relaxation 
was first quantified by Chandrasekhar in 1942 [54|. Two-body relaxation becomes important on 
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timescales longer than the local relaxation time. The evolution of a globular cluster over these 
timescales can still be described (at least formally) by the Boltzmann equation but with a collisional 
term added to the right-hand side. Equation [B] then takes on the form: 

f + .-v/-v^-| = r[/]. (11) 

Equation [TT] is sometimes called the collisional Boltzmann equation and the term r[/] describes the 
effect of two-body (and in principle higher-order) interactions on the distribution function. Practi- 
cally speaking it is not possible to evaluate r[/] analytically and various numerical approximations 
must be employed. Approaches include the Fokker-Planck method, where r[/] is approximated 
in the weak scattering limit by an expansion in powers of the phase-space parameters; the Monte 
Carlo method, where r[/] is approximated by a Monte Carlo selection of weak encounters over 
a time shorter than the relaxation time; or direct N-body integration, where rather than solving 
Equation 1111 the orbits of each star in the cluster are explicitly integrated. Each method has its 
strengths and weaknesses and will be discussed further in Section [5j 

Thermodynamically speaking, strongly self-gravitating systems have a negative heat capacity. 
This can be understood by relating the kinetic energy of the system to a dynamical temperature: 

1 3 

-mvl,^ = l^keT (12) 

|40| . This, in combination with the virial theorem, can be used to define a total internal energy 
for the cluster: 

E = -^NkBT (13) 

where N is the number of bodies in the system. Finally, this can be used to calculate a heat 
capacity: 

dB^_3 

dT 2 ^ ' 

Since all the constants on the right hand side of Equation [T3] are positive, C is always negative. 
A negative heat capacity means that heating a self-gravitating system actually causes it to lose 
energy. For a core-halo star cluster, the core is strongly self-gravitating while the halo is not, so 
the halo acts as a heat bath for the core. Any perturbation in which the core becomes dynamically 
hotter than the halo causes energy to flow into the halo. The negative heat capacity means that 
the core becomes even hotter, increasing the flow of energy to the halo in a runaway process. 
This causes the core to contract, formally to a singularity. The runaway process is called the 
gravothermal catastrophe and the consequent shrinking of the core is called core collapse. It 
affects all self-gravitating systems and was first noted in the context of star clusters by Anotnov 
[TT] . In equal- mass systems core collapse will occur after 12 — 20trh [356] but may be accelerated 
in systems with a spectrum of masses due to mass segregation. Core collapse not only appears 
in analytic, theoretical models, but has also been repeated in a variety of numerical simulations 
such as the model shown in Figure [51 [205] . Furthermore, the Harris catalogue lists several Galactic 
globular clusters that from there surface brightness profiles are thought to have experiences a core 
collapse event |151| . 

Core collapse can be halted, at least temporarily, by an energy source in the core. For stars 
(also self-gravitating systems) this energy source is nuclear burning. In star clusters tightly bound 
binaries perform a similar role. Stars in the core scatter off these binaries, gaining kinetic energy 
at the expense of the orbital energy of the binary, and the collective effect causes the core to 
return to equilibrium or even to re-expand. In analogy to nuclear burning in stars, this process in 
globular clusters is called "binary burning" . The binaries taking part in binary burning may be 
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Figure 5: Lagrange radii indicating the evolution of a Plummer model globular cluster for an 
A'^-body simulation and a Monte Carlo simulation. The radii correspond to radii containing 0.35, 
1, 3.5, 5, 7, 10, 14, 20, 30, 40, 50, 60, 70, and 80% of the total mass. Figure taken from Joshi et 
al. [205] . 
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either primordial (binaries where the stars were born bound to each other) or dynamicaUy formed 
by a variety of interactions that wih be discussed in Subsection 14.31 

Stars can escape from a star cluster if they gain a velocity greater than the cluster's escape 
velocity, {vD = —AU/M where U is the potential energy of the star cluster and M its total mass. 
Using the virial theorem it is possible to show that {vD — 4{v'^) where 2(u^) is the RMS velocity 
in the cluster [40]. There are two means through which a star can reach the escape velocity. The 
first is ejection where a single strong interaction, such as occurs during binary burning, gives the 
star a sufficient velocity impulse to exceed Ve- This process is highly stochastic. The second is 
evaporation where a star reaches escape velocity due to a large number of weak encounters during 
the relaxation process. Relaxation tends to maintain a local Maxwellian in the velocity distribution 
and, since a Maxwellian distribution always has a fraction 7 — 7.38 x 10^'^ stars with v > 2wrms, 
there are always stars in the cluster with a velocity above the escape velocity. Thus, it is the fate 
of all star clusters to evaporate. The evaporation time can be estimated as: 

tcvup ~ — = 136trlx- (15) 

7 

This is much longer than a Hubble time so few globular clusters are directly affected by evaporation. 
Evaporation can, however, be accelerated by the presence of a tidal field. In this case stars 
whose orbits extend beyond rt are stripped from the cluster and lost. Tidal dynamics are more 
complicated than a simple radial cutoff would imply and detailed prescriptions taking into account 
orbital energy and angular momentum as well as a time-varying field for star clusters in elliptical 
orbits are necessary to capture all of the important processes [3671 196] . It seems likely that the 
ultimate fate of most globular clusters is destruction due to tidal effects [1311 1119) . Since most 
compact objects are likely to live deep in the core of globular clusters where ejection will normally 
be due to violent interactions, the details of tidal stripping are unlikely to be critical for the 
treatment of relativistic binaries in star clusters. 
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3 Observations 



By their very nature, relativistic binaries are compact and faint, so that observations of these 
systems are difficult unless they are in an interacting phase. Furthermore, observations of binaries 
that have segregated to the centers of clusters can be complicated by crowding issues. Nonetheless, 
observations across a broad spectrum using a variety of ground- and space-based telescopes have 
revealed populations of binaries and their tracers in many Milky Way and extra-galactic globular 
clusters. 

Although they generally are not binaries, blue straggler stars may be used as tracers of the 
underlying binary population and the dynamical interaction rate within individual globular clus- 
ters. Blue stragglers are found on the main sequence above and to the left of the main sequence 
turn off in the CMD of a globular cluster (see Figure [5]). As their name suggests, these stars 
are too hot and too massive to still be on the main sequence if they are coeval with the stellar 
population of the cluster. They are thought to be rejuvenated through mass transfer, merger, or 
direct collision of binaries |104l [751 12991 12171 1292] . As such, they can be interpreted as tracers 
of binary interactions within globular clusters. However, many of the populations that have been 
observed exhibit a bimodal radial distribution [103] |388] [102] |329l |379l IMl [Ml ES] for which 
the inner population can be interpreted as tracing dynamical interactions and the outer population 
as representative of the primordial binary population. Recent observations of the blue straggler 
populations of 13 globular clusters indicates a correlation between the specific frequency of blue 
stragglers and the binary fraction in the globular cluster [354) . This supports observations which 
also seem to suggest that binary coalescences are the dominant formation mechanism for blue 
stragglers in globular clusters [234]. 

Blue stragglers are some of the most visible and populous evidence of the dynamical interactions 
that can also give rise to relativistic binaries. Current observations that have been revealing the 
blue straggler population and their variable counterparts (SX Phe stars) are the ACS Survey 
of Galactic Globular Clusters [3331 and the Cluster AgeS Experiment (CASE) [51]. For a good 
description of surveys that use far-ultraviolet in detecting these objects, see Knigge j216) . For 
somewhat older but still valuable reviews on the implications of blue stragglers on the dynamics 
of globular clusters, see Hut |189] and Bailyn [16]. 

The globular cluster population of white dwarfs can be used to determine the ages of globular 
clusters [272], and so they have been the focus of targeted searches despite the fact that they are 
arguably the faintest electromagnetically detectable objects in globular clusters. These searches 
have yielded large numbers of globular cluster white dwarfs. For example, a recent search of w Cen- 
tauri has revealed over 2000 white dwarfs |273) , while Hansen et al. [149] have detected 222 white 
dwarfs in M4. Deep ACS observations of NGC 6397 [326] have identified a substantial population 
of approximately 150 white dwarfs [363) . Dieball et al. [83] have found approximately 30 white 
dwarfs and about 60 probably cataclysmic variables and white dwarf binaries in M80. In general, 
however, these searches uncover single white dwarfs. Optical detection of white dwarfs in binary 
systems tends to rely on properties of the accretion process related to the binary type. Therefore, 
searches for cataclysmic variables (CVs, see Section 13. ip generally focus on low-luminosity X-ray 
sources [IMl [HHI MM and on ultraviolet-excess stars [HOI IMl UMl [HI [112], but these systems 
are usually a white dwarf accreting from a low mass star. The class of "non-flickerers" which have 
been detected recently [62] 1368] have been explained as He white dwarfs in binaries containing 
dark CO white dwarfs [92] [MT] 114811364] . 

Pulsars, although easily seen in radio, are difficult to detect when they occur in hard binaries, 
due to the Doppler shift of the pulse intervals. Thanks to an improved technique known as 
an "acceleration search" [259] . which assumes a constant acceleration of the pulsar during the 
observation period, more short orbital period binary pulsars are being discovered [45] [47] HT] 169] 
I109| I113| 1315] . For a good review and description of this technique, see Lorimer [238] . The 



14 



10 



1 

F606-F8U 



Figure 6: CMD of NGC 6205 from the Hubble Space Telescope Advanced Detector for Surveys. 
Different populations are identified as horizontal branch (HB), red giant branch (RGB), main 
sequence turn off (MSTO), and blue stragglers (BS). Figure taken from Leigh, Sills, & Knigge [235) . 
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progenitors of the ultracompact millisecond pulsars (MSPs) are thought to pass through a low- 
mass X-ray binary phase [TT] [Hni [Ml ISISl 13^ (LMXBs, see Section E^I ■ These systems are very 
bright and when they are in an active state, they can be seen anywhere in the Galactic globular 
cluster system. There are, however, several additional LMXBs that are currently quiescent |1401 
I164| 13751 1145| (144] (qLMXBs). As these systems turn on, they can be added to the list of known 
LMXBs, which is currently at 15 [162) . Additional evidence of a binary spin-up phase for MSPs 
comes from measurements of their masses, which indicate a substantial mass-transfer phase during 
the spin-up. Several observed globular cluster MSPs in binary systems are seen to have masses 
above the canonical mass of 1.4 [111] , 

Although there are many theoretical predictions of the existence of black holes in globular 
clusters (see, e.g., |265|I307|I264| [70]). there are very few observational hints of them. Measurements 
of the kinematics of the cores of M15 |116|I143] . NGC 6752 ^9], and w Centauri [2_82' provide some 
suggestions of a large, compact mass. One proposed explanation would be a large black hole of 
~ lO'^M© in the cluster core. A black holes in this mass range are much larger than stellar mass 
black holes (10 — IOOMq) but much smaller than the supermassive black holes (> 10^) found 
in galactic centers. As such they are called intermediate mass black holes (IMBH) and may 
potentially form through mergers of stars or stellar-mass black holes in globular clusters. However, 
these observations can also be explained without requiring an IMBH [2481 1290] and the existence 
of such objects is questionable. VLA observations of M80, M62, and M15 do not indicate any 
significant evidence of radio emission, which can be used to place some limits on the likelihood of 
an accreting black hole in these clusters. However, uncertainties in the expected gas density makes 
it difficult to place any upper limits on a black hole mass }25j . The unusual millisecond pulsar 
in the outskirts of NGC 6752 has also been argued to be the result of a dynamical interaction 
with a possible binary intermediate mass black hole in the core [60] . If the velocity dispersion in 
globular clusters follows the same correlation to black hole mass as in galactic bulges, then there 
may be black holes with masses in the range 1 — lO"^Af0 in many globular clusters [391] . Recent 
Hubble Space Telescope observations of the stellar dynamics in the core of 47 Tuc have placed an 
upper bound of 1500 Af© for any intermediate mass black hole in this cluster j253l . Stellar mass 
black hole binaries may also be visible as low luminosity X-ray sources, but if they are formed in 
exchange interactions, they will have very low duty cycles and hence are unlikely to be seen |2Q6]. 
For a good recent review on neutron stars and black holes in globular clusters, see Rasio et al. 321]. 
Heinke |162] has an excellent review of X-ray observations of Galactic globular clusters, including 
CVs, LMXBs and qLMXBs, MSPs, and other active binaries. 

Recent observations and catalogs of known binaries are presented in the following Subsections. 

3.1 Cataclysmic variables 

Cataclysmic variables (CVs) are white dwarfs accreting matter from a companion that is usually a 
dwarf star or another white dwarf. They have been detected in globular clusters through identifica- 
tion of the white dwarf itself or through evidence of the accretion process. White dwarfs managed 
to avoid detection until observations with the Hubble Space Telescope revealed photometric se- 
quences in several globular clusters |63l IMl EHU |323l [3241 |325l [3681 UMI ■ Spectral identification of 
white dwarfs in globular clusters has begun both from the ground with the VLT [2711 1272] and in 
space with the Hubble Space Telescope [62l [92l [3681 1273] . With spectral identification, it will be 
possible to identify those white dwarfs in hard binaries through Doppler shifts in the H/3 line. This 
approach has promise for detecting a large number of the expected double white dwarf binaries 
in globular clusters. Photometry has also begun to reveal orbital periods [276] 191] 1207] of CVs in 
globular clusters. 

Accretion onto the white dwarf may eventually lead to a dwarf nova outburst. Identifications of 
globular cluster CVs have been made through such outbursts in the cores of M5 [252] , 47 Tuc [288] , 
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NGC 6624 [Ml], M15 [342], and M22 [9l|4T]. With the exception of VlOl in M5 [252], original 
searches for dwarf novae performed with ground based telescopes proved unsuccessful. This is 
primarily due to the fact that crowding obscured potential dwarf novae up to several core radii 
outside the center of the cluster [3391 1341) . Since binaries tend to settle into the core, it is not 
surprising that none were found significantly outside of the core. Subsequent searches using the 
improved resolution of the Hubble Space Telescope eventually revealed a few dwarf novae close 
to the cores of selected globular clusters ^3381 ElOl [3441 [342| [9]. To date, there have been 14 
found [2971 1337] . using the Hubble Space Telescope and Las Campanas Observatory (CASE). 

A more productive approach has been to look for direct evidence of the accretion around the 
white dwarf. This can be in the form of excess UV emission and strong Ha emission [99] 11391 
I218| I219| l80] from the accretion disk. This technique has resulted in the discovery of candidate 
CVs in 47 Tuc EH EE], M92 [M], NGC 2808 80 , NGC 6397 [62l|92l[368], and NGC 6712 [10(1] . 
The accretion disk can also be discerned by very soft X-ray emissions. These low luminosity X- 
ray binaries are characterized by a luminosity Lx < lO'^'*'^ erg/s, which distinguishes them from 
the low-mass X-ray binaries with Lx > 10'^® erg/s. Initial explanations of these objects focused 
on accreting white dwarfs [19] , and a significant fraction of them are probably CVs [375} 1381] . 
There have been 10 identified candidate CVs in 6752 [303], 15 in 6397 [58], 19 in 6440 [302], 2 
in a; Gen [HT], 5 in Terzan 5 [TSD], 22 in 47 Tuc [HO], 5 in M80 [M], 7 in M54 [SH], 2-5 in 
NGC 288 [HI], 4 in M30 [MI], 4 in NGC 2808 [355], 1 in M71 [HI], and 1 in M4 [25]. However, 
some of the more energetic sources may be LMXBs in quiescence [375] . or even candidate QSO 
sources [26l. 

The state of the field at this time is one of rapid change as Chandra results come in and optical 
counterparts are found for the new X-ray sources. A living catalog of CVs was created by Downes 
et al. [S3], but it has been allowed to lapse and is now archival |85[ 186] . It may still be the best 
source for confirmed CVs in globular clusters up to 2006. 

3.2 Low- mass X-ray binaries 

The X-ray luminosities of low-mass X-ray binaries are in the range Lx ^ 10^^ - 10'^® erg/s. The 
upper limit is close to the Eddington limit for accretion onto a neutron star, so these systems must 
contain an accreting neutron star or black hole. All of the LMXBs in globular clusters contain an 
accreting neutron star as they also exhibit X-ray bursts, indicating thermonuclear flashes on the 
surface of the neutron star [203] . Compared with ~ 100 such systems in the galaxy, there are 15 
LMXBs known in globular clusters. The globular cluster system contains roughly 0.1% of the mass 
of the galaxy and roughly 10% of the LMXBs. Thus, LMXBs are substantially over-represented 
in globular clusters. 

Because these systems are so bright in X-rays, the globular cluster population is completely 
known - we expect new LMXBs to be discovered in the globular cluster system only as quiescent 
(qLMXBs) and transient systems become active. The 15 sources are in 12 separate clusters. Five 
have orbital periods greater than a few hours, six ultracompact systems have measured orbital 
periods less than 1 hour, and four have undetermined orbital periods. A member of the ultracom- 
pact group, 4U 1820-30 (X1820-303) in the globular cluster NGC 6624, has an orbital period of 11 
minutes [359] . This is one of the shortest known orbital periods of any binary and most certainly 
indicates a degenerate companion. The orbital period. X-ray luminosity, and host globular clusters 
for these systems are given in Tabled) 

The improved resolution of Chandra allows for the possibility of identifying optical counter- 
parts to LMXBs. If an optical counterpart can be found, a number of additional properties and 
constraints for these objects can be determined through observations in other wavelengths. In par- 
ticular, the orbital parameters and the nature of the secondary can be determined. So far, optical 
counterparts have been found for X0512-401 in NGC 1851 [ED], X1745-203 in NGC 6440 [377] . 
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X1746-370 in NGC 6441 X1830-303 in NGC 6624 ^Mi, X1832-330 in NGC 6652 [77l[T59]. 
X1850-087 in NGC 6712 [Ml [HI IMI], X1745-248 in Terzan 5 [US], and both LMXBs in NGC 
7078 [131 1385] . Continued X-ray observations will also further elucidate the nature of these sys- 
tems [275] . 

Table 1: Low-mass X-ray binaries in globular clusters: Host clusters and LMXB properties. 



LMXB Name 


Cluster 


Lx 


Porh 


Ref. 






(xlO^e erg/s) 


(hr) 




X0512-401 


NGC 1851 


1.9 


0.28 


1771 13481 13921 


X1724-307^ 


Terzan 2 


4.3 




[771 1348] 


XI 730-335 


Liller 1 


2.2 




[771 1348] 


X1732-304 


Terzan 1 


0.5 




[7711348] 


X1745-203 


NGC 6440 


0.9 


8.7 


[771 1348117] 


CX-2 


NGC 6440 


1.5-0.4 


0.96 


[resiis] 


XI 745-248 


Terzan 5 






[77] 


J17480-2446 


Terzan 5 


37 


21.25 


[351(253] 


X1746-370 


NGC 6441 


7.6 


5.70 


[771 13011 1348] 


X1747-313 


Terzan 6 


3.4 


12.36 


[771130111348] 


X1820-303 


NGC 6624 


40.6 


0.19 


[771130111348] 


X1832-330 


NGC 6652 


2.2 


0.73 


[771 1301] 


X1850-087 


NGC 6712 


0.8 


0.33 


[771 13011 1348] 


X2127-hll9-l 


NGC 7078 


3.5 


17.10 


[771130111348] 


X2127-M19-2 


NGC 7078 


1.4 


0.38 


[77] 1301] 13851 182] 



" Sidoli et al. [348] give X1724-304. 

The 15 bright LMXBs are thought to be active members of a larger population of lower luminos- 
ity quiescent low mass X-ray binaries (qLMXBs) [386^. Early searches performed with ROSAT data 
(which had a detection limit of 10'^^ erg/s) revealed roughly 30 sources in 19 globular clusters [203] . 
A more recent census of the ROSAT low luminosity X-ray sources, published by Verbunt [374] , 
lists 26 such sources that are probably related to globular clusters. Recent observations with the 
improved angular resolution of Chandra have begun to uncover numerous low luminosity X-ray 
candidates for CVs HIl IISSl IISII EMI MB M M EIZl l3n31l3S2]- For a reasonably complete 
discussion of recent observations of qLMXBs in globular clusters, see Verbunt and Lewin [375' or 
Webb and Barret [381] and references therein. Tabled] lists the 27 qLMXBs known to date. 

3.3 Millisecond pulsars 

The population of known millisecond pulsars (MSPs) is one of the fastest growing populations of 
relativistic binaries in globular clusters. Several ongoing searches continue to reveal millisecond 
pulsars in a number of globular clusters. Previous searches have used deep multi-frequency imaging 
to estimate the population of pulsars in globular clusters [113J . In this approach, the expected 
number of pulsars beaming toward the earth, A'puis, is determined by the total radio luminosity 
observed when the radio beam width is comparable in diameter to the core of the cluster. If the 
minimum pulsar luminosity is L^i^ and the total luminosity observed is itot, then, with simple 
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Table 2: Quiescent Low-mass X-ray binaries in globular clusters: Host clusters, ID numbers, and 
references. 



qLMXB Name 


Cluster 


Ref. 


171411-293159 


NGC 6304 




11451 


171421-292917 


NGC 6304 




1145) 


171433-292747 


NGC 6304 




11451 


No. 3 


LO Cen 


|145| 


|161j 


Ga 


M13 


11451 


|161) 


X5 


47 Tuc 


11451 


|161) 


X7 


47 Tuc 


Il45| 


|161) 


No.24 


M28 


Il45| 


|161| 


A-1 


M30 


|145| 


|161] 


U24 


NGC 6397 


11451 


|161J 


CX2 


MSG 


|145| 


|161j 


CX6 


MSG 


|145| 


|161j 


C2 


NGC 2808 




|145j 


16 


NGC 3201 




|145| 


CXI 


NGC 6440 




[161] 


CX2 


NGC 6440 




|161j 


CX3 


NGC 6440 




|161j 


CX5 


NGC 6440 




[161] 


CX7 


NGC 6440 




|161| 


CXIO 


NGC 6440 




|161j 


CX12 


NGC 6440 




|161j 


CX13 


NGC 6440 




|161j 


W2 


Terzan 5 




|161| 


W3 


Terzan 5 




[161] 


W4 


Terzan 5 




|161] 


W8 


Terzan 5 




|161j 


J180916-255623 


NGC 6553 




[144] 
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assumptions on the neutron star luminosity function, 



-tJmin nl (,-tJtot/ -t^min ; 

In their observations of 7 globular clusters, Fruchter and Goss have recovered previously known 
pulsars in NGC 6440, NGC 6539, NGC 6624, and 47 Tuc [HS]- Their estimates based on Equa- 
tion (|16p give evidence of a population of between 60 and 200 previously unknown pulsars in 
Terzan 5, and about 15 each in Liller 1 and NGC 6544 [113j . Additional Fermi LAT work uses 
7-ray emission to estimate total numbers of ~ 2600 — 4700 MSPs in the Galactic globular cluster 
system, which corresponds to estimates from encounter rates [Bl- 

Current searches include the following: Arecibo, which is searching over 22 globular clus- 
ters |171) : Green Bank Telescope (GBT), which is working both alone and in conjunction with 
Arecibo [202) 1171] : the Giant Metrewave Radio Telescope (GMRT), which is searching over about 
10 globular clusters |110) : and Parkes, which is searching over 60 globular clusters [55]. Although 
these searches have been quite successful, they are still subject to certain selection effects such as 
distance, dispersion measure, and acceleration in compact binaries |46| . For an excellent review 
of the properties of all pulsars in globular clusters, see the review by Camilo and Rasio [46] and 
references therein. The properties of known pulsars in binary systems with orbital period less than 
one day are listed in Table [S] which has been extracted from the online catalog maintained by 
Freire [213] • 

With the ongoing searches, it can be reasonably expected that the number of millisecond pulsars 
in binary systems in globular clusters will continue to grow in the coming years. 



3.4 Black holes 

There have been no confirmed observations of black hole binaries in Galactic globular clusters. 
All of the Galactic globular cluster high luminosity LMXBs exhibit the X-ray variability that is 
indicative of nuclear burning on the surface of a neutron star. It is possible that some of the 
recently discovered low luminosity LMXBs may house black holes instead of neutron stars |375j , 
it is more likely that they are simply unusual neutron star LMXBs in quiescence |386J . Finally, 
there is very circumstantial evidence for the possible existence of an intermediate mass black hole 
(IMBH) binary in NGC 6752 based upon an analysis of the MSP binary PSR A 59, 60 . Further 
evidence for black holes in globular clusters will probably have to wait for the next generation of 
gravitational wave detectors. 



3.5 Extragalactic Globular Clusters 

Recently the X-ray observatories XMM- Newton, Chandra, and Swift have been combined with HST 
observations to search for extragalactic globular cluster binaries. These have yielded a number of 
systems found in M31, M104, NGC 1399, NGC 3379, NGC 4278, NGC 4472, and NGC 4697. 

Three super soft X-ray sources have been found in M31 globular clusters [169ll298|ll70j . These 
are assumed to be from classical novae, and two have been associated with optical novae [1691 1170] . 
Stiele et al. [360] have found 36 LMXBs and 17 candidate LMXBs that are co-located with M31 
globular clusters, while Peacock et al. [291] claim 41 LMXBs in 11% of the M31 globular clusters. 
Barnard et al. [24] have found 4 black hole x-ray binaries associated with M31 globulars. Finally, 
there is the suspected intermediate mass black hole in Gl [lO] [78]. 

Searches for bright X-ray sources in other nearby galaxies have yielded LMXB's thought to 
contain black holes based on their luminosity. There are 2 in NGC 1399 [3461 11961 1247] . and 
Paolillo et al. [287] estimate that 65% of globular clusters in NGC 1399 host LMXBs. Brassington 
et al. [43] have found 3 LMXBs in globular clusters associated with NGC 3379, of which 1 is a 
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Table 3: Short orbital period binary millisecond pulsars in globular clusters. Host clusters and 
orbital properties. See http : //www . naic . edu/~pf reire/ GCpsr . html for references to each pulsar. 



Pulsar 


spin 


Cluster 


Porh 


e 


A/2 




(ms) 




(days) 




(Me) 


J0024-7204I 


3.485 


47 Tuc 


0.229 


< 0.0004 


0.015 


J0023-7203J 


2.101 


47 Tuc 


0.121 


< 0.00004 


0.024 


J0024-7204O 


2.643 


47 Tuc 


0.136 


< 0.00016 


0.025 


J0024-7204P 


3.643 


47 Tuc 


0.147 


— 


0.02 


J0024-7204R 


3.480 


47 Tuc 


0.066 


— 


0.030 


J0024-7203U 


4.343 


47 Tuc 


0.429 


0.000015 


0.14 


J0024-7204V 


4.810 


47 Tuc 


0.227 


— 


0.34(?) 


J0024-7204W 


2.352 


47 Tuc 


0.133 


— 


0.14 


J0024-7204Y 


2.196 


47 Tuc 


0.522 


— 


0.16 


J1518+0204C 


2.484 


M5 


0.087 


— 


0.038 


J1641+3627D 


3.118 


M13 


0.591 


— 


0.18 


J1641+3627E 


2.487 


M13 


0.118 


— 


0.02 


J1701-3006B 


3.594 


M62 


0.145 


< 0.00007 


0.14 


J1701-3006C 


3.806 


M62 


0.215 


< 0.00006 


0.08 


J1701-3006E 


3.234 


M62 


0.16 


— 


0.035 


J1701-3006F 


2.295 


M62 


0.20 


— 


0.02 


B1718-19 


1004.03 


NGC 6342 


0.258 


< 0.005 


0.13 


J1748-2446A 


11.563 


Terzan 5 


0.076 


— 


0.10 


J1748-2446M 


3.569 


Terzan 5 


0.443 


— 


0.16 


J1748-2446N 


8.667 


Terzan 5 


0.386 


0.000045 


0.56 


J1748-24460 


1.677 


Terzan 5 


0.259 


— 


0.04 


J1748-2446P 


1.729 


Terzan 5 


0.363 


— 


0.44 


J1748-2446V 


2.073 


Terzan 5 


0.504 


— 


0.14 


J1748-2446ae 


3.659 


Terzan 5 


0.171 


— 


0.019 


J1748-2446ai 


21.228 


Terzan 5 


0.851 


0.44 


0.57 


J1748-2021D 


13.496 


NGC6440 


0.286 


— 


0.14 


J1807-2459A 


3.059 


NGC 6544 


0.071 


— 


0.010 


J1824-2452G 


5.909 


M28 


0.105 


— 


0.011 


J1824-2452H 


4.629 


M28 


0.435 


— 


0.2 


J1824-2452I 


3.932 


M28 


0.459 




0.2 


J1824-2452J 


4.039 


M28 


0.097 




0.015 


J1824-2452L 


4.100 


M28 


0.226 




0.022 


J1836-2354A 


3.354 


M22 


0.203 




0.020 


J1905+0154A 


3.193 


NGC6749 


0.813 




0.09 


J1911-5958A 


3.266 


NGC 6752 


0.837 


< 0.00001 


0.22 


J1911+0102A 


3.619 


NGC 6760 


0.141 


< 0.00013 


0.02 


J1953+1846A 


4.888 


M71 


0.177 




0.032 


B2127+11C 


30.529 


M 15 


0.335 


0.681 


1.13 


J2140-2310A 


11.019 


M30 


0.174 


< 0.00012 


0.11 



presumed black hole. Fabbiano et al. [97] have found 4 LMXBs in globular clusters associated with 
NGC 4278, all of which are bright enough to be black holes. Maccarone and collaborators [24511246] 
have found 2 LMXBs containing black holes in NGC 4472 globular clusters. Kim et al. [210] claim 
75 globular clusters LMXBs in NGC 3379, 4278, and 4697. Finally, Li et al. [236] have found 41 
X-ray sources that are presumably LMXBs in globular clusters in M104. 
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4 Relativistic Binaries 



Relativistic binaries are binary systems containing two degenerate or collapsed objects and an 
orbital period such that they will be brought into contact within a Hubble time. (Note that 
this definition also includes binaries which are already in contact.) Frequently, in the process 
of becoming a relativistic binary, a binary will exist with a single degenerate or collapse object 
and a normal star. These systems are tracers of relativistic binaries. Outside of dense stellar 
clusters, most relativistic binary systems arise from primordial binary systems whose evolution 
drives them to tight, ultracompact orbits. Dynamical processes in globular clusters can drive 
wide binary systems toward short orbital periods and can also insert degenerate or collapsed stars 
into relativistic orbits with other stars. Before addressing specific evolutionary scenarios, we will 
present the generic features of binary evolution that lead to the formation of relativistic binaries. 

4.1 Binary evolution 

The evolution of a binary system of two main-sequence stars can significantly affect the evolution 
of both component stars if the orbital separation is sufficiently small. If the orbital period is less 
than about 10 days, tidal interactions will have circularized the orbit during the pre- and early 
main-sequence phase [1321 13891 1390) . Both stars start in the main sequence with the mass of the 
primary Mp and the mass of the secondary Ms, defined such that Mp > Mg. The binary system 
is described by the orbital separation r, and the mass ratio of the components q = Mg/Mp. The 
gravitational potential of the binary system is described by the Roche model where each star 
dominates the gravitational potential inside regions called Roche lobes. The two Roche lobes meet 
at the inner Lagrange point along the line joining the two stars. Figure [7] shows equipotential 
surfaces in the orbital plane for a binary with q — 0.4. If either star fills its Roche lobe, matter 
will stream from the Roche lobe filling star through the inner Lagrange point to the other star in a 
process known as Roche lobe overflow (RLOF). This mass transfer affects both the evolution of the 
components of the binary as well as the binary properties such as orbital period and eccentricity. 

Roche lobe overflow can be triggered by the evolution of the binary properties or by evolution 
of the component stars. On the one hand, the orbital separation of the binary can change so that 
the Roche lobe can shrink to within the surface of one of the stars. On the other hand, stellar 
evolution may eventually cause one of the stars to expand and fill its Roche lobe. When both stars 
in the binary are main-sequence stars, the latter process is more common. Since the more massive 
star will evolve first, it will be the first to expand and fill its Roche lobe. At this stage, the mass 
exchange can be conservative (no mass is lost from the binary) or non-conservative (mass is lost). 
Depending on the details of the mass exchange and the evolutionary stage of the mass-losing star 
there are several outcomes that will lead to the formation of a relativistic binary. The primary star 
can lose its envelope, revealing its degenerate core as either a helium, carbon-oxygen, or oxygen- 
neon white dwarf; it can explode as a supernova, leaving behind a neutron star or a black hole; 
or it can simply lose mass to the secondary so that they change roles. Barring disruption of the 
binary, its evolution will then continue. In most outcomes, the secondary is now the more massive 
of the two stars and it may evolve off the main sequence to fill its Roche lobe. The secondary can 
then initiate mass transfer or mass loss with the result that the secondary can also become a white 
dwarf, neutron star, or black hole. 

The relativistic binaries that result from this process fall into a number of observable categories. 
A WD-MS or WD- WD binary may eventually become a cataclysmic variable once the white dwarf 
begins to accrete material from its companion. If the companion is a main-sequence star, RLOF 
can be triggered by the evolution of the companion. If the companion is another white dwarf, 
then RLOF is triggered by the gradual shrinking of the orbit through the emission of gravitational 
radiation. Some WD-WD cataclysmic variables are also known as AM CVn stars. If the total mass 
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of the WD-WD binary is above the Chandrasekhar mass, the system may be a double degenerate 
progenitor to a Type la supernova. 

The orbit of a NS-MS or NS-WD binary will shrink due to the emission of gravitational 
radiation. At the onset of RLOF, the binary will become either a low-mass X-ray binary (if the 
donor star is a white dwarf or main sequence star with M < 2Mq), or a high-mass X-ray binary 
(if the donor is a more massive main-sequence star). These objects may further evolve to become 
millisecond pulsars if the neutron star is spun up during the X-ray binary phase [74 ll320j . A NS-NS 
binary will remain virtually invisible unless one of the neutron stars is observable as a pulsar. A 
BH-MS or BH-WD binary may also become a low- or high-mass X-ray binary. If the neutron star 
is observable as a pulsar, a BH-NS binary will appear as a binary pulsar. BH-BH binaries will be 
invisible unless they accrete matter from the interstellar medium. A comprehensive table of close 
binary types that can be observed in electromagnetic radiation can be found in Hilditch |172j . 
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Figure 7: Cross section of equipotential surfaces in the orbital plane of a binary with q = 0.4. The 
values of the potential surfaces are 5.0, 3.9075, 3.8, 3.559, 3.2, 3.0, and 2.8. The units have been 
normalized to the orbital separation, so a = 1. 

The type of binary that emerges depends upon the orbital separation and the masses of the 
component stars. During the evolution of a 10 Mq star, the radius will slowly increase by a factor 
of about two as the star progresses from zero age main sequence to terminal age main sequence. 
The radius will then increase by about another factor of 50 as the star transitions to the red giant 
phase, and an additional factor of 10 during the transition to the red supergiant phase. These last 
two increases in size occur very quickly compared with the slow increase during the main-sequence 
evolution. Depending upon the orbital separation, the onset of RLOF can occur any time during 
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the evolution of the star. Mass transfer can be divided into three cases related to the timing of the 
onset of RLOF: 



Case A: If the orbital separation is small enough (usually a few days) , the star can fill its Roche lobe 
during its slow expansion through the main-sequence phase while still burning hydrogen 
in its core. 

Case B: If the orbital period is less than about 100 days, but longer than a few days, the star will 
fill its Roche lobe during the rapid expansion to a red giant with a helium core. If the 
helium core ignites during this phase and the transfer is interrupted, the mass transfer is 
case BB. 

Case C: If the orbital period is above 100 days, the star can evolve to the red supergiant phase 
before it fills its Roche lobe. In this case, the star may have a CO or ONe core. 

The typical evolution of the radius for a low metallicity star is shown in Figure HI Case |X] mass 
transfer occurs during the slow growth. Case IbI during the first rapid expansion, and Case ICl during 
the final expansion phase. The nature of the remnant depends upon the state of the primary during 
the onset of RLOF and the orbital properties of the resultant binary depend upon the details of 
the mass transfer. 
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Figure 8: Evolution of the radius for a 10 M0 star with a metallicity oi Z — 0.001. Figure taken 
from Pfahl et al. [295] . 
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4.2 Mass transfer 



Although there are still many unanswered theoretical questions about the nature of the mass 
transfer phase, the basic properties of the evolution of a binary due to mass transfer can easily 
be described. The rate at which a star can adjust to changes in its mass is governed by three 
time scales. The dynamical time scale results from the adiabatic response of the star to restore 
hydrostatic equilibrium, and can be approximated by the free fall time across the radius of the 
star, 
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where M and R are the mass and radius of the star. The thermal equilibrium of the star is restored 
over a longer period given by the thermal time scale 
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where L is the luminosity of the star. Finally, the main-sequence lifetime of the star itself provides 
a third time scale, which is also known as the nuclear time scale: 
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The rate of mass transfer/loss from the Roche lobe filling star is governed by how the star's 
radius changes in response to changes in its mass. Hjellming and Webbink [1781 describe these 
changes and the response of the Roche lobe to mass changes in the binary using the radius-mass 
exponents, C, = In i?/d In M, for each of the three processes described in Equations (|171 1181 fT9|) 
and defining 

CL = (l + g)^ (20) 
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for the Roche lobe radius-mass exponent. If Cl > Cdyn, tbe star cannot adjust to the Roche lobe, 
then the mass transfer takes place on a dynamical time scale and is limited only by the rate at 
which material can stream through the inner Lagrange point. If ^dyn > Cl > Cth, then the mass 
transfer rate is governed by the slow expansion of the star as it relaxes toward thermal equilibrium, 
and it occurs on a thermal time scale. If both Cdyn and Cth are greater than Cl, then the mass loss 
is driven either by stellar evolution processes or by the gradual shrinkage of the orbit due to the 
emission of gravitational radiation. The time scale for both of these processes is comparable to the 
nuclear time scale. A good analysis of mass transfer in cataclysmic variables can be found in King 
et al. [113 . 

Conservative mass transfer occurs when there is no mass loss from the system, and therefore 
all mass lost from one star is accreted by the other star. During conservative mass transfer, the 
orbital elements of the binary can change. Consider a system with total mass M = Mi -t- M2 and 
semi-major axis a. The total orbital angular momentum 
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is a constant, and we can write a cx {M1M2) 
values by a subscript i, we find: 
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Using Kepler's third law and denoting the initial 



MuM2^ 



1 3 



M1M2 



(22) 



26 



Differentiating Equation (|22p and noting that conservative mass transfer requires Mi — 



— -M2 gives: 



P _ 3Mi (Ml - M2) 
P " M1M2 



(23) 



Note that if the more massive star loses mass, then the orbital period decreases and the orbit 
shrinks. If the less massive star is the donor, then the orbit expands. Usually, the initial phase of 
RLOF takes place as the more massive star evolves. As a consequence, the orbit of the binary will 
shrink, driving the binary to a more compact orbit. 

In non-conservative mass transfer, both mass and angular momentum can be removed from the 
system. There are two basic non-conservative processes that are important for the formation of 
relativistic binaries - the common-envelope process and the supernova explosion of one component 
of the binary. The result of the first process is often a short-period, circularized binary containing 
a white dwarf. Although the most common outcome of the second process is the disruption of the 
binary, occasionally this process will result in an eccentric binary containing a neutron star or a 
black hole. 

Common envelope scenarios result when one component of the binary expands so rapidly that 
the mass transfer is unstable and the companion becomes engulfed by the donor star. This can 
happen if the mass transfer rate is so great that it exceeds the Eddington mass accretion rate of the 
accretor, or when the donor expands past the outer Lagrange point [28411383] . The companion can 
then mechanically eject the envelope of the donor star. There are two approaches to determining 
the outcome of the process of ejection of the envelope. 

The most commonly used approach is the "a-prescription" of Webbink |382) . in which the 
energy required to eject the envelope comes from the orbital energy of the binary and thus the 
orbit shrinks. The efficiency of this process determines the final orbital period after the common 
envelope phase. This is described by the efficiency parameter 



where A_Ebind is the binding energy of the mass stripped from the envelope and Ai^orb is the 
change in the orbital energy of the binary. The result of the process is the exposed degenerate 
core of the donor star in a tight, circular orbit with the companion. This process can result in a 
double degenerate binary if the process is repeated twice or if the companion has already evolved 
to a white dwarf through some other process. If enough orbital energy is lost it can also lead to a 
merger of the binary components. A brief description of the process is outlined by Webbink [382], 
and a discussion of the factors involved in determining acE is presented in Sandquist et al. |332] . 

The other approach that is used is the "7-prescription" of Nelemans et al. [278] . In this 
approach, other energy sources such as tidal heating or luminosity of the donor may assist in the 
unbinding of the envelope. The material lost by the envelope is assumed to carry away angular 
momentum and reduce the total angular momentum of the system such that: 



where Ji and Jf are the initial and final angular momenta of the binary, M^ot is the total mass of 
the binary prior to mass loss, and AM is the mass lost in the ejection of the envelope. 

The effect on a binary of mass loss due to a supernova can be quite drastic. Following Pad- 
manabhan ^286j . this process is outlined using the example of a binary in a circular orbit with 
a semi-major axis a. Let v be the velocity of one component of the binary relative to the other 
component. The initial energy of the binary is given by 
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Following the supernova explosion of Mi, the expanding mass shell will quickly cross the orbit of 
Af2, decreasing the gravitational force acting on the secondary. The new energy of the binary is 
then 

1 / MnsMs \ 2 GMnsM2 



= 2 Imn^TI^ j ^' - 

where Mns is the mass of the remnant neutron star. We have assumed here that the passage 
of the mass shell by the secondary has negligible effect on its velocity (a safe assumption, see 
Pfahl et al. |295) for a discussion), and that the primary has received no kick from the supernova 
(not necessarily a safe assumption, but see Davies and Hansen [7T or Pfahl et al. |296] for an 
application to globular cluster binaries). Since we have assumed that the instantaneous velocities 
of both components have not been affected, we can replace them by = G (Mi + M2) /a, and so 



2a \Mns + M-^ 

Note that the final energy will be positive and the binary will be disrupted if A/ns < (1/2) (Mi + 
M2). This condition occurs when the mass ejected from the system is greater than half of the 
initial total mass, 

AM > i (Ml + M2) , (29) 

where AAf = Mi — M^^. If the binary is not disrupted, the new orbit becomes eccentric and 
expands to a new semi-major axis given by 
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and orbital period 



Note that we have not included any mention of the expected velocity kick that newly born neutron 
star or black hole will receive due to asymmetries in the supernova explosion. These kicks can be 
quite substantial, up to several hundred kilometers per second and, at least for observed pulsars, 
seem to be drawn from a Maxwellian distribution with a peak at 265 km s~^ [179] . In most cases, 
the kick will further increase the likelihood that the binary will become unbound, but occasionally 
the kick velocity will be favorably oriented and the binary will remain intact. If the kick is higher 
than the escape velocity of the cluster, it will also remove the remnant from the system. This 
mechanism may be very important in depleting the numbers of neutron stars and black holes in 
globular clusters. 

We have seen that conservative mass transfer can result in a tighter binary if the more massive 
star is the donor. Non-conservative mass transfer can also drive the components of a binary 
together during a common envelope phase when mass and angular momentum are lost from the 
system. Direct mass loss through a supernova explosion can also alter the properties of a binary, 
but this process generally drives the system toward larger orbital separation and can disrupt the 
binary entirely. With this exception, the important result of all of these processes is the generation 
of tight binaries with at least one degenerate object. 

The processes discussed so far apply to the generation of relativistic binaries anywhere. They 
occur whenever the orbital separation of a progenitor binary is sufficiently small to allow for 
mass transfer or common envelope evolution. Population distributions for relativistic binaries are 
derived from an initial mass function, a distribution in mass ratios, and a distribution in binary 
separations. These initial distributions are then fed into models for binary evolution such as 
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Star Track |32] or SeBa |309[ I279j in order to determine rates of production of relativistic binaries. 
The evolution of the binary is often determined by the apphcation of some simple operational 
formulae such as those described by Tout et al. [370] or Hurley et al. [183) . For example, Hils, 
Bender, and Webbink [177] estimated a population of close white dwarf binaries in the disk of 
the galaxy using a Salpeter mass function, a mass ratio distribution strongly peaked at 1, and a 
separation distribution that was flat in ln(a). Other estimates of relativistic binaries differ mostly 
by using different distributions [3T | [T95 l [279 ] [277] . 

4.3 Globular cluster processes 

In the galactic field stellar densities are low enough that stars and binaries rarely encounter each 
other. In this environment binaries evolve in isolation with their properties and fate determined by 
solely by their initial conditions and the processes described in Sections 14. II and [4.21 In star cluster 
stellar densities are much higher and encounters between stars, stars and binaries, and binaries 
and binaries become increasingly common. Such encounters can affect a binary's parameters and 
dramatically alter the evolution that would be expected if the same binary were isolated. Some 
encounter outcomes that are of particular interest for relativistic binaries include: 

• Reduction or increase of the period during distant encounters. 

• Exchange interactions where binary membership is changed in close encounters. 

• Binary formation during strong few-body interactions. 

• Encounter-induced binary mergers. 

Different types of interactions can produced each of these outcomes and we will describe these 
briefly in the terms normally used by stellar dynamicists - the number of bodies and type of 
objects involved in the interaction. In all cases these interactions must be distinguished from the 
very weak interactions that drive two-body relaxation of clusters as described in Section 12.31 The 
interactions affecting binary parameters result from close encounters, the outcomes of which cannot 
be described statistically using the language of relaxation theory. Note also that the term binding 
energy is normally spoken of as though it were a positive quantity. Thus binaries with the largest 
(most negative) binding energy are the most tightly bound. 

4.3.1 Single-Single Interactions 

As the name suggests, single-single interactions are encounters between two single stars. If the 
periastron of the encounter is sufficiently small, the stars may excite tidal oscillations in each other 
at the cost of some of their relative kinetic energies. If sufficient kinetic energy is dissipated, the 
stars can become bound and form a new binary. 

The exact nature of the oscillations excited in the stars are not important, only that they 
dissipate sufficient kinetic energy. Furthermore, to obtain a bound orbit it is only necessary to 
dissipate sufficient kinetic energy that the total energy at apocenter is reduced to less than zero. 
The basic condition for tidal capture is: 

AE^i + AEt2 > ^t^vl (32) 

where A_Eti,T2 is the energy associated with the tidal oscillation in each star, fi is the reduced 
mass, and Va is the relative velocity at apocenter (or infinity) [HE]. Only a fraction ~ (va/vp)^, 
where Vp is the velocity at pericenter, of the energy must be converted to tidal oscillations in order 
for a capture to occur. For example, an encounter with a pericenter velocity of ^ lOOkms"^ in 
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a cluster with a velocity dispersion of ~ lOkms^^ only needs to dissipate ~ 1% of the apocenter 
kinetic energy in order for a capture to occur. This mechanism favours the creation of binaries in 
encounters with large ratios between apo and pericenter velocity. Thus it tends to produce highly 
eccentric binaries with very small pericenter separations. Such binaries would then be expected to 
circularize over time due to further tidal effects. 

This process was once thought to be the dominant path for dynamically creating binaries in star 
clusters [98l HO] since two-body interactions are more likely than higher-order encounters [356] . It 
has been realized, however, that this is not the case since stars must approach very closely, to within 
a few stellar radii, in order for tidal oscillations to dissipate enough energy for a capture to occur 
[98l I233| 1311) . In these situations it is actually more likely, depending upon the stellar equation of 
state, for the stars to merge rather that form a binary |3221 1255[ \T9\ 11941 1320j . Furthermore, due 
to the difhculty in exciting tides in degenerate matter, tidal capture is unlikely to form binaries 
from two compact objects. 

4.3.2 Three-body Interactions 

When three single-stars interact it is possible for one of the stars to gain kinetic energy and, if the 
relative kinetic energy of the remaining two stars is sufficiently reduced, they can form a bound pair. 
In general the least massive body will gain the highest velocity, since it is the easiest to accelerate, 
and the two more massive bodies will remain (this can also be seen as a consequence of energy 
equipartition) . Heggie |153) showed that if Eb is the binding energy of the new binary, the rate of 
three-body binary formation is oc i^j^ . Thus three-body interactions tend to favour the creation 
of binaries with low binding energies. The binding energy can be increased by later encounters (see 
Section [4.3. 3 1) and there is about a 10% probability that a binary formed by three-body interactions 
will gain sufficient binding energy to survive permanently [15^ 1356] . 

Three-body binary formation interactions are less common than tidal interactions since the 
probability of a close encounter between three stars is smaller than the probability of an encounter 
between two. Thus stellar densities need to be higher for three-body binary formation to be efficient 
than for tidal interactions |134l 11881 1356] . However, because it is not necessary to excite tides, 
pericenter passages need not be as close in order to form binaries during these interactions and 
mergers are a less common outcome of three-body interactions. It is worth noting that, in equal- 
mass systems, neither tidal nor three-body interactions are likely to produce many binaries over 
the lifetime of a star cluster [1881 11341 l40l 1356] . This is not necessarily the case for clusters with 
a mass function due to the very high central densities and strong encounters that can be achieved 
by mass segregation. Nevertheless, primordial binaries and the interactions involving them are 
critical for producing large numbers of relativistic binaries in globular clusters. 

4.3.3 Binary - Single Interactions 

Binary-single interactions, although formally still three-body interactions, differ because two of 
the stars are already bound. In such an encounter several outcomes are possible depending on the 
relative kinetic energy of the star and binary and on the binding energy of the binary itself. 

If the kinetic energy of the star is less than the binding energy of the binary, energy equipartition 
requires that the star be accelerated and the binding energy of the binary be increased (the period 
will become shorter). If, however, the kinetic energy of the star is greater than the binding energy 
of the binary, the star will donate kinetic energy to the binary and the binding energy will decrease 
(the period will become longer). For equal-mass systems this introduces a simple yet important 
classification for binaries. If \Eb\ < ksT (Section 12. 3|) then the binary is said to be "soft" and, on 
average, will tend to have its period increased by encounters (softening). If \Ei,\ > ksT than the 
binary is said to be "hard" and will tend to have its period decreased by encounters (hardening). A 
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hard binary gains on average (AEt) ^ 0.2 — QAEb per interaction assuming a Maxwellian velocity 
distribution [1531 1192] (although this is only fully valid in the case of very hard binaries |188) ). 
Since the encounter rate is proportional to the semi-major axis of the binary (oc !/£'(,), the energy 
of a hard binary increases by AE^^rix ~ O-QkbT/t^ix per relaxation time [40]. This leads to the 
"Heggie- Hills Law" |153[ 1174] that states "hard binaries get harder while soft binaries get softer" . 
The end result of binary hardening can be a merger while the end result of binary softening can 
be the disruption of the binary. 

In a multi-mass system the division between hard and soft binaries is not so clear since the 
relative energies depend upon the specific masses of all the bodies involved. For an individual 
encounter between a binary with member masses mi and m2 and a single star of mass ms travelling 
at velocity w, it is possible to define a critical velocity |190] : 

2 mi + 1712 + ms 

Vc = 7 r (33) 

TO3(toi + m2) 

such that for Vc < v the binary cannot capture the single star and the star can disrupt the binary 
(softening) and for Vc > v then the binary can capture the single star and cannot be disrupted 
in the process (hardening). Thus the distinction between hard and soft binaries still exists in 
multi-mass systems. It has been shown by Hills that it is also possible to use the ratio of orbital 
speeds to define whether the binary will gain or lose energy |173 \ 

The reference to capture alludes to another important process that can occur in binary-single 
interactions: exchange. In an exchange interaction one of the original binary members is replaced 
by the single star so the binary membership changes. As with other encounters, equipartition 
of energy favours the ejection of the lowest-mass object and exchange encounters are a way of 
introducing massive objects into binaries. It has been shown that in the limit of m^ ^ mi or 
m2, the probability for a massive object to be exchanged into the binary is ~ 1 '176'. Thus it is 
possible to create a compact binary from a binary where neither member was originally compact 
entirely by dynamical means. This may be particularly effective for BHs since they are the most 
massive objects in evolved star clusters [350] . 

4.3.4 Binary - Binary Interactions 

There are many possible outcomes for binary-binary interactions, particularly in the multi-mass 
case, that depend rather sensitively upon the initial conditions of the encounter. Therefore a 
quantitative theory of these interactions is lacking. For sufficiently distant encounters both binaries 
will "see" each other as single centres of mass and a hardening or softening interaction will occur in 
both binaries, depending on their binding energies and relative kinetic energy. There are also some 
general results for close encounters where all stars have equal mass. If one of the binaries is much 
harder than the other, it can be exchanged into the other binary as a single star. This produces 
a hierarchical three-body system, the future of which will be decided by further interactions [356] . 
One or both binaries can also be disrupted during the encounter. Numerical experiments show 
that at least one of the binaries is disrupted in '^88% of cases |260|I261|I262II115] . It is also possible 
for one or, more rarely, both binaries to exchange members with each other. Thus binary-binary 
interactions can provide all of the effects of binary-single interactions but can result in extensive 
binary destruction as well. 

Higher-order interactions (those involving more than four bodies) are also possible but will be 
quite rare and are even less amenable to quantitative description than are binary-binary encoun- 
ters. 
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5 Dynamical Evolution 



The evolution of the compact binary population in a globular cluster depends on the interplay 
between the global evolution of the cluster as described in Section 12.31 the few-body interactions 
described in Section 14.31 and the binary stellar evolution described in Section 14.11 Simulations 
that hope to describe the compact binary population of a globular cluster must be able to couple 
these processes and describe their mutual interaction. 

The evolution of globular clusters is primarily governed by point-mass gravitational attraction 
between individual stars and thus their overall structure can be described in terms of classical 
Newtonian A^-body dynamics. Other processes, however, also play a role. Stellar evolution ob- 
viously affects the number and properties of compact objects and compact binaries so describing 
it accurately is important for producing the correct initial compact population. Stellar evolution 
proceeds on the same timescales as relaxation and core-collapse in globular clusters and can result 
in significant mass-loss from individual stars. This mass-loss can change the mass of the cluster 
enough to affect its dynamical evolution. Consequently the global dynamical evolution of a glob- 
ular cluster is coupled to nuclear burning in individual stars and this must be taken into account. 
Furthermore, stellar evolution affects the orbital parameters of the binary population and, as we 
have seen in Section [2.31 binaries can be an energy source of for the cluster and can slow down or 
halt core collapse. Thus the details of both binary evolution and few-body interactions can affect 
the global dynamics of a globular cluster. Finally, the loop is closed by the fact that the global 
dynamics a globular cluster can affect the properties of the binary population by bringing stars 
and binaries close enough that they can interact and perturb each other. We are interested in the 
binary population that is the end result of this complex interplay of different processes operat- 
ing on very different spacial and temporal scales. Simulating star cluster dynamics is a problem 
that has a long history and has promoted contact between various branches of astrophysics. The 
MODEST (Modeling DEnse STellar systems) collaboration, a collection of various working groups 
from different backgrounds, maintains a website that provides some up-to date information about 
efforts to model the dynamical and stellar evolution of star clusters |270l. 

5.1 Star Cluster Simulation Methods 

Broadly speaking, there are three approaches that can be used to simulate the dynamical evolution 
of star clusters. The first, direct iV-body integration, treats the dynamical evolution of a star 
cluster by numerically solving Newton's equations. The second, the distribution function method, 
represents the the star cluster using a 6A^ dimensional distribution function and then describes its 
evolution using a collisional form of the Boltzmann equation. Both of these methods have their 
strengths and weaknesses. Several authors have compared these techniques (e.g. |154| 1212| 1158| 
1371) ) and in general find that they agree well. The summary of the MODEST-2 meeting |353j gives 
a fairly recent review of efforts to implement these techniques and use them in simulations. The 
third method, the encounter rate technique, does not self-consistently simulate the evolution of star 
clusters but can be uses a static cluster model combined with individual scattering experiments to 
estimate the effect of dense stellar systems on a specific population. In the next three sections we 
outline these techniques. We then conclude the chapter by summarizing some of the most recent 
results about compact binary populations in star clusters derived using these methods. 

5.1.1 Direct A^-body Integration 

Direct A^-body simulations are the most conceptually straightforward method of simulating star 
clusters. Each star is explicitly represented by a massive (non-test) particle and the gravitational 
interactions are calculated by numerically integrating the 3A^ coupled differential equations of 
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motion in classical Newtonian gravity: 




\fj - fi|3 



(34) 



Here rrij is the mass of particle j and fij are the positions of particles i and j respectively. A 
great deal of general information on direct A'^-body methods can be found in the review article by 
Spurzem |357j and Aarseth's book on A^-body simulations f?. 

There are currently two major families of codes used for direct A^-body simulations: the 
Kira/STARLAB/AMUSE environment and the NBODYA series of codes. Kira integrates the 
A^-body equations of motion using a 4*''-order Hermite predictor-corrector scheme [250) . Kira is 
available as part of the star cluster evolution code STARLAB [191 , a code that also includes a 
stellar and binary evolution package, SeBa, based on the analytic tracks of Hurley et al. |183lll86j . 
and various other routines for generating initial conditions and dealing with tidal processes. Kira 
is also available in the more recent star cluster simulation environment AMUSE (Astrophysical 
Multipurpose Software Environment) |305] . AMUSE is designed to be highly modular and allows 
the user a choice between various integrators, stellar evolution routines, and prescriptions for deal- 
ing with stellar collisions. The NBODYA series of codes are a long-standing feature of the stellar 
dynamics community and have been in constant development by Aarseth since the late 1960s [1]. 
These codes employ the same basic 4*''-order Hermite predictor-corrector scheme as Kira and in- 
clude (in some versions) stellar evolution based on the same Hurley et al. analytic tracks. The 
NBODYX codes, however, treat binary evolution rather differently and are highly optimized for 
computational efficiency. Aarseth has produced two excellent reviews of the general properties 
and development of the NBODYA codes [3j|2] and includes further details in his book on A-body 
simulations [4] . A fairly complete summary of common A^-body codes and related applications can 
also be found on the NEMO website [280]. 

The primary advantage of direct A^-body simulations is the small number of simplifying as- 
sumptions made in treating the dynamical interactions. All gravitational interactions are explicitly 
integrated so the direct A'^-body method can, in principle, resolve all of the microphysics involved 
in the dynamical evolution of globular clusters. In this sense direct A^-body simulations represent 
star clusters "as they really are" . Furthermore, since all trajectories involved in any interaction are 
known, it is easy to study each interaction in detail. Finally, because direct A^-body simulations 
produce a star-by-star models of globular clusters, it is relatively straightforward to couple them 
with additional stellar physics such as stellar evolution, stellar equations of state, stellar collision 
models and so on. Therefore, within the limits of numerical error, the results of direct A^-body 
simulations are very reliable and they represent the gold standard in globular cluster modelling. 

Direct A^-body simulations are, however, very computationally expensive. Over a relaxation 
time, the computational cost of a direct A^-body simulation intrinsically scales as A^'^"^ [155] . 
Two factors of A^ come from calculating the pairwise interactions, another factor comes from the 
increasing number of interactions per crossing time as A^ increases, and a final factor of N/logN 
come from the increasing length of ij-ix compared to the crossing timescale (see Equation [3]) . 
Another problem is the very large range of scales present in globular clusters [18911155] . Distances, 
for instance, can range from a few tens of kilometers (neutron star binaries) to several tens of 
parsecs (the size of the cluster itself). The timescales can have an even larger range, from several 
milliseconds (millisecond pulsar binaries) to several Gigayears (the relaxation times of very large 
globular clusters). In the most extreme cases it may be necessary to know positions to 1 part 
in 10^* while a cluster may need to be integrated forward on timescales 10^° times shorter than 
its lifetime. Furthermore, A^-body systems are known to be chaotic 2661 1267] 1133] . This means 
that the integration schemes must be very accurate (4*''-order or better) and, consequently, they 
too are very computationally expensive. Obviously for large N these simulations rapidly become 
numerically intractable! 
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There are several methods that can be used to mitigate these problems. The Hermite scheme, 
for instance, requires only the first and second derivatives of acceleration, even though it is a 4*'*- 
order method. Since dynamical evolution proceeds more quickly in regions of high stellar density, 
the stars in cluster halos do not need to have their phase space parameters updated as frequently 
as those in cluster cores. Both STARLAB and NBODYX include individual adaptive timesteps 
to ensure that during each force calculation, only the necessary particles have their parameters 
updated. Some of the NBODYX codes also include a neighbour scheme so that not all pairwise 
forces need to be explicitly calculated at each timestep. Binaries tight enough that their internal 
evolution is largely unaffected by external forces from the cluster can also be removed from the 
global dynamical evolution and replaced with a single centre-of-mass particle. The evolution of the 
binary itself can then be followed using an efficient special-purpose code. STARLAB calculates 
the evolution of binaries using a semi-analytic Keplarian two-body code that admits perturbations 
from the cluster environment. The NBODYX codes use a regularization technique to achieve 
the same effect with the additional advantage of removing the singularity created by very close 
binaries. Both of these methods avoid the global timestep of the star cluster being reduced to that 
of its tightest binary. 

Attempts have been made to create parallel direct A^-body codes with NB0DY6-1--I- [3571 being 
the most successful example. A major problem with this type of parallization is that each particle 
in an A^-body calculation formally requires information about every other particle for each force 
calculation. This means that parallel A^-body codes are communication dominated and do not scale 
well with the number of processors. Another approach has been to use hardware acceleration. The 
GRAPE (GRavity PipE), invented by Makino [251^, is a special-purpose computer designed to 
rapidly calculate the force necessary for direct A^-body simulations. The current generation, 
the GRAPE-6, was reported to have a peak speed of 64 T-flops in 2002 \MW . A single GRAPE 6 
board, the GRAPE-6A, is also available as PCI card. The next generation GRAPE, the GRAPE- 
DR, is predicted to have a peak speed of 2 P-flops |249j . It is not yet clear when (or if) the 
GRAPE-DR will become commercially available. Another very promising possibility is hardware 
acceleration using GPUs (Graphical Processing Units) ,358. . GPUs provide a similar service to the 
GRAPE while having the advantages of being available off-the-shelf and commercially supported. 
AMUSE, NB0DY6 and NB0DY6-K-H have been modified to make use of the GPU. 

Despite developments in both hardware and software, current direct A^-body simulations are 
practically relegated to A^ < 10^, about an order of magnitude lower than the number of stars 
in a typical globular cluster. Furthermore, such simulations can take several months to run, 
making large parameter studies impossible. For this reason more approximate methods employing 
distribution functions are still widely used in stellar dynamics. 

5.1.2 Distribution Fiinction Methods 

Some of the computational problems associated with direct A'^-body systems can be avoided by 
describing globular clusters in terms of distribution functions rather than as ensembles of particles. 
The distribution function for stars of mass m at time i, fmix, v, t), is GA^-dimensional in position- 
velocity phase space. The number of stars in the phase space volume (PxcPv is given by dN = 
fmd^xd^v. This description requires that either the phase space volume d^xd^v is small enough 
to be infinitesimal yet large enough to be statistically meaningful or, more commonly, that fm be 
interpreted as the probability of finding a star of mass m in the volume d'^x d^v at time t. The time 
evolution of the cluster is modelled by calculating the time evolution of /,„ . The effect of gravity 
is divided into two components. The first is a smoothed potential, 0, that describes the effect of 
the overall gravitational field of the cluster. (t>{x) is found by solving the Poisson equation, which 
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can be written in terms of the distribution function as: 

(35) 



The second component is a coUision term, r[/], that can be though of as the effect of random close 
encounters between individual stars. In this picture (j) represents the smooth background gravita- 
tional field of the cluster while r[/] represents the granularity induced by the mass concentration 
in individual stars. The time evolution of fm is calculated using a modified form of the Boltzmann 
Equation: 

^ + v-yfrn-y^-^ = m- (36) 

In stellar dynamics Equation 131] is usually called the Fokker-Planck equation. 

Equation 1361 can be directly numerically solved if a specific expression for r[/] can be found. 
One common choice is to consider a cluster in the weak scattering limit where distant two-body 
interactions drive dynamical evolution. This is essentially equivalent to assuming the cluster evolves 
only by two-body relaxation as described in Section [273] In this case r[/] can be described by an 
expansion in the phase space coordinates (x, v G w) and takes on the form [40) : 

= - E tIt [/(-)(A-.)] + ^ E jrS^ mw){Aw.Aw,)] (37) 
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where (AX) is the expectation value of the phase space variable over some time At. In practise 
it is not efficient to solve Equations 1361 and 1371 in Cartesian coordinates, however a careful choice 
of At can alleviate this problem. The key observation is that there are two relevant timescales 
for evolution in the Fokker-Planck approximation: the crossing time, icr, that governs changes in 
position, and the relaxation time, irix, that governs changes in energy. In the case where icr ^ ^rix 
changes in position are essentially periodic and at any given time the orbit of a star can be written 
in terms of its energy, E, and the magnitude of its angular momentum vector, J. By making a 
careful choice of At and suitable assumptions about the symmetry of the potential and the velocity 
dispersion (57] 1356] it is possible to describe the evolution of the phase space structure of the cluster 
solely in terms of the orbital parameters and Equations 1361 and 1371 can be written in terms of E and 
J alone. In this representation Equation 1361 is called the orbit-averaged Fokker-Planck equation. 

Numerically solving the orbit-averaged Fokker-Planck equation is much less computationally 
intensive than direct A^-body simulations and thus can be used to model full-sized globular clusters. 
Since these simulations are much faster (taking hours or days rather than months) they are also 
useful for performing parameter space studies. It is also possible, in principle, to improve the 
accuracy of the Fokker-Planck method by including higher-order terms in Equation l37] [334 . The 
Fokker-Planck method can be generalized to include velocity anisotropics, tidal stripping, non- 
spherical potentials and cluster rotation [95] 12111 ]212J . 

The Fokker-Planck method also has some serious limitations. Each stellar mass in the cluster 
requires its own distribution function in Equation 1361 For more than a few distribution functions 
the Fokker-Planck method becomes numerically intractable so the ability of the method to resolve 
a continuous spectrum of masses limited. Since the Fokker-Planck method uses a distribution 
function, it does not automatically provide a star-by-star representation of a globular cluster. 
This means that it can be difficult to compare the results of a Fokker-Planck simulation with 
observations. It is also difficult to include stellar evolution in Fokker-Planck simulations because 
there are no individual stars to evolve. For the same reason Fokker-Planck methods cannot easily 
model strong few-body interactions [115] . interactions that are essential for describing the evolution 
of the compact binary fraction in star clusters. Therefore Fokker-Planck methods are of limited 
use in studying binary populations in star clusters. 
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Another approach to solving Equation[36]that avoids some of the problems of the Fokker-Planck 
approach is the Monte Carlo method. This method was first developed by Henon |167l 11661 1168] 
and significantly improved by Stodolkiewicz |361l 1362] . It does not use the distribution function 
explicitly but rather represents the cluster as a collection of particles just as would a direct iV- 
body simulation. Each particle is characterized by a mass, a distance from the cluster centre, and a 
tangential and radial velocity. The underlying treatment of relaxation is the same as in the Fokker- 
Planck approach but is calculated on a particle-by-particle basis. At the beginning of a Monte 
Carlo timestep the code calculates the global potential based on the current mass and radius of each 
particle. This potential is then used to generate a plane-rosette orbit for each particle, defined by 
E and J (the orbit-averaged approximation) . Next each particle has its orbit perturbed by a weak 
encounter, the parameters of which are calculated between it an a neighbouring particle. The key 
to the method is that assuming weak scattering, this perturbation is statistically representative of 
all weak encounters the particle will experience over a time At. Thus the results of the perturbation 
can be multiplied by an appropriately chosen factor to represent the total perturbation the particle 
receives over At. Once the new orbits have been found, new positions can be randomly chosen 
from a time-weighted average along the orbital path, a new potential calculated, and the process 
can start again. 

The Monte Carlo method is mathematically identical to the Fokker-Planck method but, since 
each star is individually represented, it is much easier to couple with additional physics. Each star 
in the simulation can have a different mass with no loss of efficiency so the method can resolve a 
continuous mass spectrum. Stellar evolution can be included in exactly the same way as a direct 
iV-body code. It is also relatively easy to include strong few-body interactions. For each star or 
binary at each timestep a cross-section for a strong interaction with another star or binary can 
be calculated and the encounter resolved using either analytic prescriptions or explicit few-body 
integration. Thus the Monte Carlo method has many of the advantages of the direct A^-body 
method while being much faster. Since there are a fixed number of operations per particle, the 
Monte Carlo method scales almost directly with A^ rather than the A^'^~^ dependence of a direct 
A^-body code. Thus a simulation of A^ ~ 10® can be completed in a matter of days, rather than 
the months required by the direct A^-body method. Monte Carlo codes have been shown to re- 
produce observations |156[ 11251 1126] and the results of direct A^-body codes |122l 11061 [55] quite 
well, although some discrepancies remain [371] . 

There are several Monte Carlo codes currently available, in particular that of Freitag [112] . 
that of the Northwestern Group [2051 M [Ml ES], and that of Giersz [1201 [IHl \l22\ The 
Freitag code is focused on the simulation of star clusters with a central massive black hole but 
the Northwestern and Giersz codes are both focused on general cluster dynamics and are roughly 
comparable in capabilities. 

The Monte Carlo method does have some shortcomings. It is limited by the weak scattering 
approximation and is not always reliable in the core region where interactions can be very strong. 
Furthermore, the method can only be as good as the assumptions that go into the treatment of 
few-body interactions. In order to solve the gravitational potential at each timestep efficiently, 
it is necessary to assume that the cluster is spherically symmetric, also a significant limitation. 
Nonetheless, Monte Carlo methods have proved to be quite effective for efficiently modelling star 
clusters. 

Another option for solving equation [36] is to use an analogy between between a globular cluster 
and a self-gravitating gaseous sphere [2391 1127] . These models are solved by taking moments of 
the Boltzmann equation and, like the Fokker-Planck method, can have their accuracy increased 
by including higher-order moments [334] . These are still continuum models and have difficulty 
dealing with strong few-body interactions and stellar evolution. A solution to this is the so-called 
"hybrid" code where single stars are treated by the gaseous moment equations while strong few- 
body encounters are treated using a Monte Carlo method [1281 1129] . So far this method have not 
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been used to examine compact binary populations in star clusters. 
5.1.3 Encounter Rate Techniques 

Encounter rate techniques are another method that has proven very useful for understanding the 
effect of dense stellar environments on binary populations. These methods do not deal with the 
global evolution of globular clusters but rather deal with the interactions between an environment 
and a particular stellar species. Normally a stellar density distribution is established for a globular 
cluster (or its core) and then the number of interactions per unit time is calculated for a specific 
object using cross-sections. The density distribution, interactions and outcomes can be calculated 
either analytically or numerically, depending on the complexity of the system in question. These 
methods are limited but are very useful for estimating which processes will be the most significant 
for specific stellar populations. There is no "general method" for encounter rate techniques but we 
will discuss some specific implementations in section [5.21 

5.2 Results from Models of Globular Clusters 

Over the past two decades there has been extensive work modeling the dynamical evolution of 
globular clusters. Published studies have employed encounter rate techniques in static cluster 
backgrounds, Fokker-Planck and Monte Carlo models and direct N-body simulations (e.g. [35l 

[Ml EMI [Ml EQll EIIl 13201 EMI M 

[T971 nm nm\ M\ MQ\ [Mil [MSl [THU [TM])- These models can, if properly interpreted, provide a 
wealth of information about stellar populations in globular clusters. These models, however, have 
limitations that must be understood. The dynamical approximations employed coupled with an 
imperfect understanding of stellar and binary evolution make it impossible to compare the results 
of these simulations to observations. The problem is exacerbated by the fact that many authors 
are more concerned with the development of numerical methods than on producing astrophysical 
results. Furthermore, few of these studies focus on (and some do not even model) the compact 
binary population. Those that do often focus on the effect that the compact population will have 
on the cluster dynamics rather than the effect the cluster dynamics will have on the compact binary 
population. Nevertheless many published simulations provide useful insights into the behaviour of 
compact binary populations in dense stellar systems. Here we report on these results for binaries 
where the primary is a WD, NS or BH, trying to give both a historical overview as well as details 
of the latest results. When reading past papers, attention must be paid to terminology. A neutron 
star binary, for instance, can imply a binary composed of two neutron stars or simply a binary 
where a neutron star is paired with any type of secondary depending on the authors preference. 
Here we attempt to clearly indicate the nature of the secondary in order to resolve any ambiguities. 

5.2.1 Binaries with White Dwarf Primaries 

Two kinds of binaries where the white dwarf is the most massive degenerate member have attracted 
particular attention from modelers: cataclysmic variables (CVs) and double white dwarf binaries 
(DWDs). As described in Chapter 13.11 CVs consist of white dwarfs accreting matter from a 
companion, normally either a dwarf star or another white dwarf. CVs with a dwarf star companion 
are not, strictly speaking, relativistic binaries as defined in the introduction. Unlike many double- 
degenerate binaries, however, CVs can be observed in the electromagnetic spectrum and they are 
very useful for understanding the behaviour of compact objects in globular clusters. CVs with a 
helium white dwarf companion are know as AM CVn systems and are a type of DWD. DWDs, 
which as the name suggests are binaries of two white dwarfs, are relativistic binaries but not all 
are AM CVns. The merger of DWDs can, however, be observed as a SNe la and thus they too 
have attracted some attention from globular cluster modelers. 
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Cluster 


Mtot [Mq] 


Po [Mq/pc] 


■Tc [pc] 


n [pc] 


trh [Gyr] 


w Cen 


3.34 X 10" 


4.0 X 10^ 


3.4 


74.0 


16.0 


47 Tuc 


1.06 X 10" 


5.1 X 10* 


0.59 


69.0 


8.7 



Table 4: Structural parameters for 47 Tuc and uj Cen taken from Davies and Benz [73]. Mtot is 
the total mass, po the central density, the core radius, rj the tidal radius and trh the current 
estimated half-mass relaxation time. 

White dwarfs are the lowest-mass compact stellar remnants and are not necessarily significantly 
more massive than the main sequence turn-off of present-day globular clusters. Because of this 
they may not experience as significant mass segregation as higher-mass stellar remnants and may 
not be so clearly preferred for exchange interactions. As such, modeling is crucial to understanding 
their behaviour in globular clusters. 

Lee (1987) [232] is an example of an early numerical study of the compact binary population 
in globular clusters, focusing on the effect of the binaries on the cluster dynamics. It is based on 
Fokker-Planck models using two different stellar mass components to represent main sequence stars 
and compact objects. It demonstrated that compact objects have a strong effect on the cluster 
core, implying they experience many interactions and highlighting the importance of the cluster 
environment for compact binaries. 

One of the early, specific studies of CVs in globular clusters is found in Di Stefano and Rap- 
paport (1994) [T^ who performed a large series of experiments using an encounter-rate technique. 
Specifically, they performed numerical two-body scattering experiments in static cluster back- 
grounds coupled with binary stellar evolution prescriptions in order to investigate the tidal capture 
and subsequent evolution of CVs. The cluster parameters were chosen to match the current states 
of 47 Tuc and w Cen. Since these clusters are frequently compared in studies of compact binaries, 
some of their properties are given in Table |4] (more can be found in the Harris Catalogue of globular 
clusters |151] ). The primary differences between these clusters are that w Cen is significantly more 
massive that 47 Tuc while 47 Tuc is much more concentrated and dynamically older than oj Cen. 

Di Stefano and Rappaport report 150 — 200 CVs in 47 Tuc, ^ 45 of which have an accretion 
luminosity of > 10^'^ erg s~^, and are thus visible, at the current epoch. For u Cen they predict 
100 — 150 and ^ 20 respectively. They claim a factor of < 10 enhancement in the number of CVs 
over the field population. The authors point out both that this number of CVs was barely within 
the upper limits on the number of CVs in clusters known observationally at the time and that the 
number could be even higher when binary-single and binary-binary interactions were taken into 
account. Observations of CVs in globular clusters were in their infancy at the time. However the 
over-production of CVs would remain a problem in later works. 

The next studies to touch on CVs in globular clusters were a series of papers by Davies and 
collaborators in 1995-1997 [73] [71] [72] concerned with exotic binaries produced by dynamical 
interactions. These were also encounter rate models but included binary-single interactions as well 
as tidal capture. Davies and collaborators [73[[7T] reported many single-single interactions between 
white dwarfs and main sequence stars but found that the majority of strong encounters will lead 
to a white dwarf surrounded by the envelope of a disrupted star rather than a CV. Overall Davies 
predict more CVs in 47 Tuc (smaller, concentrated, dynamically older) than than in lo Cen (larger, 
less concentrated, dynamically younger) but more merged systems than CVs in both. Davies 
cautions, however, that merging may be over-estimated during tidal capture. They also predict a 
reasonable number of DWD mergers in both clusters, most of which are above a Chandrasekhar 
mass and will thus lead to SNe la but they spend little time on this population. 

In 1997 Davies [72] produced another encounter rate study concentrating on the effect a cluster 
environment would have on a population of binaries that would become CVs due to binary stellar 
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evolution in isolation. This study observed that disruption of binaries in a dense stellar environment 
is important and can reduce the number of CVs compared to the field. Davies found that most CV 
progenitors in cluster cores with number densities > 10'^ are ionized during dynamical interactions. 
Thus, if dynamical interactions form no new CVs, the cores of galactic globular clusters will be 
depleted in CVs. CV progenitors in the low-density halos of globular clusters may survive to 
become CVs if the relaxation time of the cluster is long enough that they do not mass-segregate 
to the core. Davies proposed that comparing the CV population in the cores and halos of globular 
clusters could be used as a probe of their dynamical properties. A cluster with a core depleted 
in CVs compared to the halo will not have experienced much dynamical binary formation while a 
cluster with a large number of CVs in the core must have experienced dynamical binary formation 
at some point in its history. 

The first detailed study of DWDs in dense stellar environments became available with the 
direct N-Body simulations of Shara and Hurley (2002) |343] . They were specifically interested 
in the galactic rate of SNe la and if globular clusters could significantly enhance it. They used 
the NB0DY4 code to perform four 22 000 particle simulations with 10% primordial binaries, two 
different metallicities and detailed stellar evolution prescriptions. These simulations are in the 
open rather than globular cluster regime and were terminated after only 4.5 Gyr when they had 
lost more than 75% of their initial mass. They were run on 16 chip GRAPE-6 boards and took 
approximately five days each to complete. 

Shara and Hurley found a factor of 2-3 enhancement in the number of DWDs in a clustered 
environment compared to a similar field population. Exchange interactions seem to be particularly 
important since of the 135 DWDs produced between all four simulations, 93 are formed by ex- 
changes. The enhancement in SNe la progenitors, DWDs with a total mass above a Chandrasekhar 
mass and a gravitational wave inspiral timescale shorter than a Hubble time, is even greater. They 
find ~4 SNe la progenitors per 2000 binaries, a factor of ~15 more than in the field. Both exchange 
interactions and hardening during fiy-bys are responsible for the enhancement. Any DWD that 
went through a common envelope phase were assumed to be circularized and they found that later 
encounters did not induce significant eccentricity in the binaries. Thus they concluded that most 
DWDs in globular clusters will be circular. They find that neither normal CVs nor AM CVns are 
significantly enhanced over the field population however those in the clusters are often the result 
of exchange interactions. Thus while star clusters will not necessarily enhance the number of AM 
CVns, their properties may be different from the AM CVns in the field. This also implies that most 
CVs in star clusters are the product of dynamics and that pure binary stellar evolution channels 
for CV and AM CVn creation are shut down in clustered environments. 

Shara and Hurley (2006) [345] continued this work using a single 100 000 simulation with 5000 
primordial binaries. This simulation took six months to run to an age of 20 Gyr using a 32-chip 
GRAPE-6 board, demonstrating the strong A^-dependence of computation time in direct N-body 
simulations. The focus of this simulation was on normal CVs, 15 of which were produced during the 
course of the simulation and three of which would be visible at the current epoch. This illustrates 
the limitation of direct N-body methods. After six months of computing time the authors had only 
15 objects to analyze. Further simulations, either for verification purposes or to explore parameter 
space, are simply too computationally expensive. By contrast the encounter rate techniques were 
able to produce many simulations over a large parameter space and produce a much larger number 
of binaries for better statistics. Nonetheless the much more accurate N-body simulations are very 
important for verifying the approximate results. Shara and Hurley found that dynamical effects 
pushed several of the progenitors to the active CV state in a shorter time that would have occurred 
in the field. Indeed one of these binaries would have merged without becoming a CV if not for 
an orbital perturbation. Furthermore the simulation produced two binaries formed by exchange 
interactions that clearly have no analogs in the field. The exchange CVs are, on average, more 
massive than CVs in the field. They confirm that many CV progenitors are destroyed by dynamical 
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interactions and that a field population of stars with similar properties would produce 17 CVs over 
20 Gyr, nine of which would be visible at the current epoch. Thus Shara and Hurley predict '-^2 — 3 
times less CVs in globular clusters than in the field. 

This simulation, although primarily focused on CVs, includes some discussion of the DWD 
population and generally confirms the picture presented in Shara and Hurley 2002. In particular 
short-period (P < 1 day, comparable to the largest period of SNe la progenitors in Shara and 
Hurley 2002) DWDs are enhanced by a factor of 5-6 over the field population. Although this is 
somewhat less than the factor of ~ 15 previously reported it is still much larger than any possible 
enhancement in CVs. Shara and Hurley propose several reasons for more efficient DWD production 
in star clusters. In a DWD both components must be sufficiently massive to become white dwarfs 
within a Hubble time while only one does in a standard CV. Thus CVs are a much more common 
outcome of pure binary stellar evolution than are DWDs. DWDs are also, on average, more massive 
than CVs and thus tend to mass-segregate further towards the core and undergo more interactions. 
Finally, since each component of a DWD is likely to be more massive than the main sequence star 
in a CV and since exchanges tend to place massive objects in binaries, a DWD is a more likely 
outcome of an exchange interaction between a main sequence star and two white dwarfs than is a 
CV. 

One of the most complete projects on the compact binary population in the recent literature is 
found in a series of papers by Ivanova and collaborators (3201 1197U201II200T 119911198] . These papers 
simulate the binary - and particularly the compact binary - population of star clusters using a 
combination of Monte Carlo simulations and well-developed encounter rate techniques coupled with 
stellar and binary evolution algorithms tailored towards compact objects |183[ 1186] l33]. Ivanova 
et al. (2005a) |197) consider the effect of cluster dynamics on the initial binary fraction of globular 
clusters and show that dynamical disruption will significantly reduce the binary fraction of for 
reasonable initial conditions. Indeed they show that a globular cluster where 100% of stars are 
initially binaries will normally have a core binary fraction of only ^10% after several Gyrs. This 
implies that globular clusters may have had high binary fractions when young and thus binary- 
single and binary-binary interactions may be more important than previously thought. They also 
demonstrate both the binary burning that supports the cluster against core collapse and a build-up 
of short-period white dwarf binaries in the core as shown in Figure [9] 

The models were improved by Ivanova et al. (2006) |200j in order to investigate the white dwarf 
binary population in globular clusters. A more sophisticated, time-evolving two-zone model was 
used as a cluster background. Scattering encounters were chosen using a Monte Carlo method and 
numerical few-body integration for single-single, binary-single and binary-binary interactions was 
used to determine their outcome. They confirm that CVs are not highly over-produced in globular 
clusters relative to the field population but find that CV production is moderately efficient in 
cluster cores. In the cores they predict an enhancement over the field population of a factor of 
2-3. As with Shara and Hurley, the globular cluster CVs have different properties than field CVs. 
In particular the globular cluster CVs are slightly brighter and have a larger X-ray-to-optical flux. 
Ivanova et al. attribute this to slightly higher-mass white dwarfs in the cluster CVs, as observed 
by Shara and Hurley. The more massive white dwarfs lead to the globular cluster CVs extracting 
more energy at the same mass transfer rate as field CVs and having slightly higher magnetic fields. 
This would account for the observed X-ray excess. 

Ivanova et al. can shed more light on the dynamical processes that form CVs than previous 
encounter rate studies since they include most of the likely formation mechanisms in the same 
simulation. Their results show that some 60-70% of CVs in the core form because of dynamical 
interactions while most of the CVs in the halo originate from primordial binaries. Overall they 
find that only ^^25% of CVs formed in the clusters would have become CVs in the field. They 
conclude that tidal capture of main sequence stars is largely unimportant for CVs but that collisions 
between red giants and main sequence stars - a situation where the envelope of the red giant can 
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Figure 9: Binary period distributions from the Monte Carlo simulation of binary fraction evolution 
in 47 Tuc. The bottom panel indicates the period distribution for binaries containing at least one 
white dwarf. is the total number of binaries and A^b,p is the number of binaries per bin. Figure 
taken from Ivanova et al. |197J . 
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stripped away leaving its core as white dwarf bound to the main sequence star - can provide 
up to 15% the total number of CVs in the clusters. Exchange interactions are the most efficient 
channel, forming ~ 32% of all CVs while collisions during binary-single or binary-binary encounters 
produce another ~ 13%. The remaining CVs are produced by primordial binaries but 20% of 
these experienced hardening encounters in fly-by interactions and would not have become CVs in 
the field within a Hubble time. Finally, Ivanova et al. claim that 60 % of CVs did not form solely by 
common envelope evolution, the most common channel for field CV production and thus highlight 
the importance of cluster dynamics for the CV population. Ultimately Ivanova et al. predict 1 
detectable CV per 1000-2000 Mq of mass in Galactic globular clusters. They specifically predict 
35-40 CVs in 47 Tuc, of the same order as the 22 observed there at the time. The discrepancy 
could be due to stochasticity, some inherent deficiencies in the modeling or because not all CVs in 
47 Tuc have been observed. Observations of globular clusters in the 2000s, particularly in the UV 
[1631 ISn 1244) have indicated a larger number of CVs than previously thought. The combination 
of fewer CVs in current models and a larger number of detected CVs has gone a long way towards 
solving the "CV problem" in globular clusters. 

Unlike Shara and Hurley, Ivanova et al. find little evidence for a strong enhancement in the 
number of DWDs. This could be '"real" in the sense that their simulations model larger systems 
than Shara and Hurcly's work and thus could be taken as evidence that while DWDs are enhanced 
in open clusters, they are not in globular clusters. It could also be a result of the different 
numerical methods. The two-zone model of Ivanova et al., for example, may not be able to 
represent mass segregation as accurately as the direct N-body models of Shara and Hurley and 
thus the effect of interactions in the dense core could be underestimated. It is also possible that 
the encounter prescriptions could be missing some details that are important for DWD creation. 
More controlled comparisons of the two methods are necessary in order to definitively state which 
prediction is more accurate. Ivanova et al. do, however, find a modest enhancement in the number 
of SNe la progenitors over the field population. They predict that this will not be important for 
understanding the SNe la from Milky Way-type galaxies since only a very small fraction of their 
masses is in globular clusters. Elliptical galaxies, however, tend to have a larger fraction of their 
mass in globular clusters and in these galaxies the contribution of globular clusters to SNe la rate 
could be important. Ivanova et al. point out that DWD mergers and accretion induced collapse 
of white dwarfs in binaries also produce neutron stars and may do so without imparting a velocity 
kick. If natal kicks are large enough to remove most core-collapse neutron stars, merged white 
dwarfs could form a significant fraction of the current neutron star population in globular clusters. 

Willems et al. 2007 [387^ performed an extended analysis of the simulations of Ivanova et 
al. 2006 in order to determine the number of eccentric DWDs that may be in the low frequency 
gravitational wave band accessible to space-based gravitational wave detectors, such as the Laser 
Interferometer Space Antenna (LISA). Taking the DWD eccentricities from the lO"^ to 13"^ Gyr 
of the Ivanova et al. simulations, Willems et al. generate the probability distribution for eccentric 
DWDs in Galactic globular clusters shown in Figure [101 Unlike Shara and Hurley, these results 
indicate that there will be a small but non-negligible population of eccentric DWDs in the Galactic 
globular cluster system (see Table [5]). They will have no analogs in the field since field DWDs are 
expected to be circularized by common envelope evolution. Thus a detection of an eccentric DWD 
by LISA is both probable and a clear indication of a dynamical origin. 

The consensus that emerges from the simulations is that dynamical destruction of potential 
CVs and dynamical formation of new CVs are equally important in globular clusters. Therefore 
the number of CVs per unit mass is unlikely to be much higher in globular clusters than it is in the 
field. Globular cluster CVs may, however, be slightly more massive than field CVs and thus have 
slightly higher X-ray fiuxes. Overall current simulations seem to be able to re-produce the observed 
number of CVs in clusters reasonably well but still tend to slightly over-predict the population. 
Whether this is due to a failing in the simulations or observational limitations is not yet clear. 
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Figure 10: Statistical properties of eccentric DWDs (e > 0.01) in the LISA band {P < 5000s) for 
three different cluster models. Hatched histograms are for the entire eccentric DWD population 
while shaded histograms are fore the DWDs with at least two harmonics individually resolvable 
by LISA with an S/N > 8 (5 year mission, 10 kpc distance). The PDFs are normalized to unity. 
From left to right: frequency, eccentricity, total mass, mass ratio and lifetime before merger. Plot 
taken from Willems et al. 2007 [387] . 
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Pz 


Typical 


0.39 


0.32 


0.06 


0.01 


47 Tuc 


1.90 


0.84 


0.55 


0.29 


Terzan 5 


1.60 


0.80 


0.48 


0.22 



Table 5: (A^) is the mean number of eccentric binaries per cluster and Pn is the Poisson probability 
that at least A^ eccentric DWDs are present at any given time. Data taken from Willems et al. 
2007 (3871. 
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It seems likely that the DWD population in globular clusters is larger per unit stellar mass than 
that in the field however it is not clear if the enhancement is by factors of ten or only factors of a 
few. The disagreement between various simulations could be due to differing numerical methods 
or because scaling the DWD population from open to globular clusters is not straightforward. 
Systematic comparison of the different simulation methods will be necessary in order to resolve 
this unambiguously. 

5.2.2 Binaries with Neutron Star Primaries 

Neutron star binaries were the first binaries with degenerate companions to be specifically inves- 
tigated in globular clusters. LMXBs, neutron stars accreting from a companion as described in 
Chapter 13.21 have attracted particular interest since results from 1975 |208l [56] gave indications 
that they are over-abundant in globular clusters. These papers immediately suggested this was a 
result of dynamical encounters and by the end of the 1970s a debate over which type of encounter, 
whether tidal capture by MS stars [HH], collisions with RG stars |365] . exchange interactions |175| 
or three-body interactions [347], was already underway. Like CVs, LMXBs are not necessarily 
relativistic by our definition but are electromagnetically active and thus useful for calibration 
purposes. The ultracompact systems discussed in Chapter 13.21 often called UCXBs, are almost 
certainly double-degenerate and may consist of a neutron star accreting from a white dwarf. 

Neutron stars are, on average, more massive than the main sequence turn-off of old globular 
clusters and also more massive than white dwarfs. As such they should both mass-segregate to 
cluster cores and be preferentially exchanged into binaries during interactions. Thus globular 
clusters may produce a large number of double neutron star binaries (DNSs), the mergers of which 
will be excellent sources of gravitational radiation for ground-based detectors and may also be a 
source of short gamma-ray bursts [142] . Neutron star production in globular clusters is, however, by 
no means certain. Neutron stars are rare and thus primordial DNSs are not common. Furthermore 
large natal kicks can both break up primordial DNSs and remove neutron stars from globular 
clusters. Also, neutron stars are only a factor of two to three times more massive than the average 
cluster mass, leaving the extent of their mass segregation and their probability of exchange unclear. 

Neutron stars are low-luminosity and, if not members of an X-ray binary, are normally detected 
as MSPs. Their presence was long-expected in globular clusters due to the large population of 
LMXBs that are thought to be MSP progenitors but were not confirmed until 1987 [243]. Since 
all MSPs are thought to need a spin-up phase during accretion in a binary system, even single 
MSPs can be used as a probe of binary evolution. As discussed in Chapter 13. 3( these are now 
known to be quite common in globular clusters. Of more interest for gravitational wave detection 
are the many binaries containing MSPs and particularly the double MSP binaries (DMSPs) that 
have been detected in some globular clusters. 

Originally tidal capture or collisions were the favoured method for producing X-ray binaries. 
When the LMXB enhancement was first discovered, globular clusters were thought to have very 
few primordial binaries and thus there would be little opportunity for exchange interactions. Since 
three-body interactions tend to produce soft binaries they were largely ignored. During the 1980s 
studies employing static cluster models and encounter rate estimates based either on analytic 
estimates (e.g. [1931 13731 [T8l 1376) ) or few-body scattering experiments (e.g. |2231 1224] ) were able 
to roughly re-produce the X-ray binary population in globular clusters using only tidal capture 
and collisions. They also indicated that hardening encounters might be important in bringing 
dynamically formed binaries into the mass-transferring, and thus X-ray, regime. Verbunt and 
Meylan |376] showed that mass segregation has a significant effect on the production of CVs and 
X-ray binaries, an early indication that globular cluster structure needs to be treated in some detail 
rather than as a homogeneous static background. 

As stellar modeling improved, the central role of tidal capture began to be called into question. 
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Several authors demonstrated that close encounters between MS or RG stars lead to mergers 
rather than bound pairs |322[ 12551 and thus will not produce X-ray binaries. It was also 
realized that, since dynamical interactions in globular cluster cores can destroy binaries [197 , a 
low binary fraction in old Milky Way globular clusters did not necessarily imply a low primordial 
binary fraction. Thus, particularly in young clusters, interactions between stars and binaries are 
likely to be important. 

One of the earliest systematic attempts to simulate double-degenerate binaries in globular clus- 
ters was the Rappaport et al. (1989) |318) study of the effect of a dense stellar environment on a 
wide binary containing an MSP with a low-mass white dwarf companion. The study utilized an 
encounter rate technique where a large number of numerical binary-single scattering experiments 
were conducted in a static cluster background with parameters chosen to represent uj Gen and 
47 Tuc. They found that binary-single interactions were capable of introducing significant eccen- 
tricity into some of the binaries, enough to effect their gravitational wave inspiral timescales. It 
also appeared that disruption interactions were unlikely, at least for binaries with P < 300 days 
so hard double-degenerate binaries should be able to survive in globular cluster cores. The results 
were marginally consistent with the available observations but could not be fully reconciled with 
predicted eccentricities as a function of age as estimated by spin-down observations. 

Sigurdsson and Phinney [3521 1351] conducted similar scattering experiments to those performed 
by Rappaport et al. |315] but with tighter binaries and a much extended mass spectrum. They 
discovered that exchanges were even more important than previously thought, and indeed could be 
the dominant type of interaction between massive stars and moderately hard binaries. A secondary 
effect of exchange is that when a light star is swapped for a massive one, the binary will have a 
larger semi-major axis for the same binding energy, which increases the cross section of the binary 
to undergo further interactions. They also found that three-body interactions with very hard 
binaries tended to produce a merger between two of the members and accretion onto the neutron 
star during such an event could plausibly create a MSP. Indeed Sigurdsson and Phinney argued 
that such interactions are sufficient to explain the entire MSP population of globular clusters. They 
also showed that a small number of low mass DMSPs can be formed in this way. 

Sigurdsson and Hernquist |349j considered the role of binary-binary interactions in forming 
DMSPs. They proposed a mechanism where two neutron star-main sequence binaries interacted 
and both main sequence stars were disrupted. The result was an DNS surrounded by an envelope 
that could be accreted to spin up both members to produce a DMSP. The study employed both 
point-mass scattering experiments and an SPH representation of the main sequence stars, when 
appropriate. The model appeared plausible but has not factored into later studies and would be a 
rather rare, if interesting, outcome of binary-binary interactions. 

Davies and collaborators [73] [71] specifically investigated the population of LMXBs in globular 
cluster in the same simulations described in Chapter 15. 2. II Davies and Benz (1995) [73 considered 
the evolution of 1000 binaries injected into various cluster cores modelled as King models with a 
concentration parameter Wq between 3 and 12 - a range that covers the concentrations of galactic 
globular clusters. In this study they found that the number of LMXBs produced depends on 
the IMF, assumptions about the how easily the envelope can be ejected during common envelope 
evolution, assumptions about stellar evolution and the concentration of the cluster. Any of these 
parameters can affect some types of binaries by factors of several. 

Davies (1995) 71 concentrated on the blue straggler and compact binary population of 47 Tuc 
and ui Gen. Particular interest was shown in comparing the efficiency of single-single and binary- 
single encounters in producing various kinds of binaries. Here the effect of different cluster models 
became clear. For the model of uj Gen (massive, low concentration) blue stragglers are produced 
significantly more efficiently by binary-single than by single-single interactions. By contrast in 
47 Tuc (lower mass, higher concentration) binary-single and single-single interactions are of com- 
parable significance in blue straggler formation. In both systems binaries containing neutron stars 
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seem to be more likely to be formed by single-single interactions with an main sequence star rather 
than binary-single interactions. Like the CVs, however, binaries so-formed, tend to lead to dis- 
ruption of the main sequence star and a "smothering" of the neutron star with the main sequence 
star's envelope. This may lead to the formation of a Thorne-Zytkow object. The neutron star 
will eventually accrete the matter which could lead to black hole production (if enough matter is 
accreted) or an MSP but will not lead to an X-ray binary. A smothered system can also be the 
result of a binary-single interactions but is a less common. For the most optimistic set of parame- 
ters Davies predicts a 4:1 ratio of smothered to non-smothered neutron star binaries in u Cen and 
a 5:1 ratio in 47 Tuc. 

Davies predicted that uj Cen would form more X-ray binaries than in 47 Tuc since the IMF in 
uj Cen model was skewed towards higher masses and thus produced a larger number of neutron 
stars. By contrast more CVs were predicted in 47 Tuc. Both clusters were predicted to produce 
more CVs than LMXBs which was difficult to reconcile with observations at the time. The number 
of predicted LMXBs was consistent with contemporary observations of 47 Tuc and uj Cen (the 
possibility of one or two in each at the current epoch). 

A new model for neutron star-white dwarf binaries was introduced by Rasio et al. (2000) 
|320) who proposed that short-period MSP-white dwarf binaries can be formed by exchanges with 
primordial binaries followed by a common envelope phase. In this model a ~l-3 M© main sequence 
star in a binary acquires an neutron star companion through an exchange interaction and forms a 
new binary at a separation of 0.1-1.0 AU. The main sequence star undergoes Roche lobe overflow 
and has its envelope stripped by the neutron star. Spin-up during the accretion phase leaves an 
MSP in a close orbit with a white dwarf. Rasio et al. use Monte Carlo scattering simulations 
to predict the MSP-white dwarf binary population produced by this scenario and compare it to 
the population of short-period MSP-white dwarf binaries observed in 47 Tuc. The results, given 
in Figure lll[ show very good agreement with the observations and confirm the importance of 
interactions for the compact binary population in globular clusters. 

The tail end of the systems in group B of Figure [11] represents the neutron star-white dwarf 
binaries that are in very short period orbits and are undergoing a slow inspiral due to gravitational 
radiation. These few binaries can be used to make an order of magnitude estimate on the population 
of such objects in the galactic globular cluster system. If we consider that there are two binaries 
with orbital period less than 2000 s out of ^ 10^ Mq in 47 Tuc, and assume that this rate is 
consistent throughout the globular cluster system as a whole, we find a total of ^ 60 such binaries. 
This is a crude estimate and may be optimistic considering the results of Ivanova et al. |200l I199j 
but is an example of how simulations can be extended to produce gravitational wave detection 
estimates. 

Ivanova et al. (2005b) [201' Re-introduced the plausibility of a collision between a red giant 
and a neutron star forming an eccentric neutron star- white dwarf UCXB. The result is based on 
hydrodynamical simulations |237j that show, while a common envelope phase is likely after such a 
collision, much of the envelope is ejected to infinity and an eccentric UCXB is a common outcome 
of the interaction. Ivanova et al. use the results from 40 SPH simulations and simple, analytic 
cluster modeling to estimate that it is possible for all UCXBs in globular clusters to be formed by 
this mechanism. 

Ivanova et al. (2008) [199] extended the results of the CV study described in Chapter 15.2.11 
to cover neutron star binaries as well. The dynamical situation for neutron star binaries slightly 
different than for CVs. About half of the neutron star binaries in the core formed as the result of 
an exchange interaction with only a few percent more formed by either tidal capture or collision. 
The rest of the binaries are primordial. Very few neutron star binaries will become X-ray binaries 
and Ivanova et al. estimate there is only a ~3% chance that any given neutron star will ever be in 
an LMXB. Furthermore, most primordial binaries that will become X-ray binaries will do so quite 
early in the life of the cluster and will have ceased to be luminous in the current epoch. Thus 
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Figure 11: Results of the Monte Carlo simulation of neutron star- white dwarf binary generation 
and evolution in 47 Tuc. Each small dot represents a binary system. The circles and error bars 
are the 10 binary pulsars in 47 Tuc with well measured orbits. Systems in A have evolved through 
mass transfer from the white dwarf to the neutron star. Systems in B have not yet evolved through 
gravitational radiation to begin RLOF from the white dwarf to the neutron star. Systems in C 
will not undergo a common envelope phase. Figure taken from Rasio et al. |320| . 
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Ivanova et al. predict that X-ray binaries observed in Galactic globular clusters today are likely 
to be products of exchange, tidal capture, or collisions. Although exchange interactions produce 
most of the dynamically formed neutron star binaries, actual X-ray binaries will be produced in 
equal numbers by all three channels. This is due to each channel having a different probability of 
producing a mass-transferring system. Finally Ivanova et al. predict ^7.5 UCXBs in the entire 
Galactic globular cluster system and ^ 180 qLMXB systems. They also make specific predictions 
that certain clusters should have at least one LMXB and all of these predictions are met. 

In addition to X-ray binaries, Ivanova et al. 2008 |199| also consider the population of MSPs 
and double neutron star binaries. Overall they have difficulty reconciling the number of X-ray 
binaries and MSPs in their simulations with the observed ratios. If they assume that all mass- 
gaining events lead to MSPs then the number produced is too high unless the detection probability 
for MSPs in clusters is only ^ 10%. They speculate that either not all accretion leads to spin up, 
that the kicks for evolution-induced collapse (EIC) supernovae are larger than they model, or their 
model for mass segregation does not work properly for these systems and more should sink to 
the core and be disrupted. They find that by adjusting their assumptions about spin-up during 
various kinds of accretion they can bring the number of MSPs down to a level that is consistent 
with observations of 47 Tuc and Terzan 5. They also note that the fraction of single MSPs to 
MSPs in binaries tends to increase with increasing cluster density. This implies that destructive 
interactions become important for MSP binaries at high density. 

They find that the production of double neutron star binaries is very inefficient. There are only 
three primordial neutron stars in the 2x10^ M© of stars that make up their entire simulation set 
while only 14 double neutron stars are formed by dynamical interactions. Given their chosen IMF 
and cluster ages, this implies that it takes ^ 500 single NSs to produce one merging double neutron 
star binary in 11 Gyrs. These results are very density dependent with the vast majority of the 
DNSs forming in the highest density simulations. They estimate that merging DNSs in globular 
clusters could provide at most 10-30% of short gamma ray bursts but note that this agrees fairly 
well with the predictions of Grindlay et al. (2006) |142j . 

This result agrees with results of Downing et al. (2010) [87] who use a similar-sized Monte 
Carlo dataset, primarily to study black hole-black hole binaries, and find no DNSs in any of their 
simulations. That DNSs are not produced in similar numbers to X-ray binaries in globular clusters 
may be a combination of the lower number of potential progenitors, kicks for both progenitors, 
and the likelihood of merging or producing a WD during interactions. 

Ivanova et al. also report a significant number of LISA sources in their simulations. They 
predict that there should be an average of 10 LISA sources per globular cluster at any given 
instant and on average 180 LISA sources are produced per globular cluster per Gyr. Scaling this 
to the Milky Way, they find that there may be up to 1500 LISA binaries, both DWD and white 
dwarf-neutron star, in the Galactic globular cluster system at the current time. They caution, 
however, that they may have overestimated the white dwarf-neutron star binary population since 
it is not fully consistent with UCXB observations. 

Banerjee and Ghosh US] have employed a novel method to investigate the compact binary 
population of globular clusters. They used the collisional Boltzmann equation as given by Spitzer 
j356' to calculate the evolution of the number distribution of binaries in a uniform, non-evolving 
background of stars. They did not employ the Fokkcr-Planck approximation but rather considered 
mean encounter rates with a randomly fiuctuating component. This is called a Wiener process and 
is a mathematical description of Brownian motion. The Wiener process means that the encounter 
rates, normally ordinary differential equations in the Boltzmann equation, now become stochastic 
differential equations and must be treated with a mathematical tool known as Ito calculus. Both 
Wiener processes and Ito calculus are described in Banerjee and Ghosh 2008b [23]. 

Banerjee and Ghosh specifically consider a bimodal population of 0.6 Mq main sequence stars 
interacting with 1.4 Mq neutron stars. They consider binary formation by tidal capture and 
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exchange, destruction by 3-body disruption of the binary or by exchange of the compact component 
for a main sequence star, and hardening by magnetic breaking and gravitational radiation. They 
are able to re-produce the build-up of short-period binaries found by other authors but, more 
interestingly they are able to compare the mean rates predicted by the Boltzmann equation to the 
effects of stochasticity in the Wiener process. Thus they are able to indicate if it is possible to 
predict the number of X-ray sources expected in a given globular cluster or if stochastic effects will 
dominate. Figure [12] shows both the mean number of X-ray binaries and the stochastic fluctuation 
expected for various cluster parameters as well as observed numbers from several Galactic globular 
clusters. Although the number of X-ray binaries in the models and observations follow the same 
trend, the stochastic fluctuations are large and indicate it is probably impossible to predict exact 
numbers of compact binaries in any given globular cluster. Thus, although the method of Banerjee 
and Gosh is in some ways less detailed than that of Ivanova et al., containing much more limited 
mass functions, a less accurate treatment of individual encounters, and much more limited stellar 
evolution, it is still valuable since it constrains the effects of stochasticity in a controllable way. 
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Figure 12: Contours of constant number of X-ray binaries as a function of the "Verbunt param- 
eters". Here F = p^r^jvc-, 7 = p/vc, p is the average stellar density, is the core radius and Vc 
the core velocity dispersion. F describes the two-body stellar encounter rate and 7 describes the 
rate of encounters between a single binary with the stellar background. The squares represent the 
location of various Galactic globular clusters in the 7 — F plane and the number in parentheses 
is the number of X-ray binaries observed in the cluster. Plot taken from Banerjee and Ghosh 
2008b [23]. 

Recently Bagchi and Ray [H] [TS] have used scattering experiments to re-investigate the eccen- 
tricities of MSP binaries in light of improved observations. They noted that observationally there 
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are three, statistically distinct, eccentricity groupings in MSP binaries, those with e '--^ 0, those 
with 2 X 10^^ < e < 0.01 and those with 0.01 < e. They then performed scattering experiments, 
utilizing the STARLAB package, to attempt to determine the origin of the different groupings. 
They found that most of the higher eccentricity binaries could be explained by exchanges or merg- 
ers during binary-single or binary-binary interactions. The intermediate eccentricities could be 
explained by fly-by encounters while the lowest eccentricities were binaries that had been circu- 
larized without experiencing fly-bys. They also found one binary that they could only explain by 
invoking resonant interactions within a hierarchical triple |222| . 

Of a more esoteric nature is the work of Trenti et al. 2008 |372| considering the number triple 
systems that might contain at least one pulsar. They studied triple formation using direct N-body 
simulations with between 512 and 19 661 bodies, both equal masses and with a mass function, and 
various initial binary fractions. They find that triple systems are about two orders of magnitude 
less common than binary systems within globular cluster cores but, given the number of MSP 
binaries currently know in globular clusters, this implies that an MSP triple should be detected 
soon. They note, however, that triples are very difficult to observationally confirm if the scales of 
the inner and outer systems are well-separated. It is also clear that if a triple with a single MSP 
will be so rare, a triple system with more than one MSP is even more unlikely. 

There are several overall conclusions that can be drawn about the population of neutron star 
binaries in globular clusters. Dynamical interactions can certainly enhance the number of X-ray 
binaries in globular clusters but the current consensus is that the enhancement will by factors of 
a few, not factors of hundreds. Recent simulations show that detailed descriptions of dynamics, 
binary stellar evolution and stellar collisions are all necessary in order for the results to be reliable. 
The same interactions that produce X-ray binaries are also effective at producing MSPs, although 
many of the resulting binaries will be subsequently disrupted and the MSPs will continue their lives 
as single objects rather than binaries. It also seems that neutron star-white dwarf binaries may be 
fairly common in globular clusters. By contrast detailed simulations show that DNS binaries are 
not common in globular clusters since few exist primordially and they are not produced efficiently 
by stellar dynamics. The individual results of these studies are highly dependent on the detailed 
treatment of various phases of binary stellar evolution as well as details of the kick distribution. 
It is worth noting that many authors have successfully predicted the number of X-ray binaries 
in globular clusters since the early 1980s. On one hand this is shows that the X-ray binary 
enhancement in globular clusters can be easily understood. On the other hand, it probably means 
that X-ray binaries are a rather generic outcome of globular cluster dynamics. Therefore observed 
numbers of X-ray binaries probably cannot tell us much about which processes are important in 
specific globular clusters. 

5.2.3 Binaries with Black Hole Primaries 

Black holes are non-luminous and can only be detected in the electromagnetic spectrum if they 
are accreting matter, normally from a binary companion. Black hole-black hole binaries (DBHs) 
can only be detected in gravitational radiation. So far there are no confirmed detections of black 
holes in globular clusters. 

Black holes are the most massive compact stellar remnants, up to 35 — SOMq |114[ 129] . and 
can, on average, be several times more massive than the average mass of stars in even a fairly 
young globular cluster. Thus black holes are strongly affected by mass segregation, may form 
Spitzer-unstable subsystems in the dense core of the cluster and will be preferentially exchanged 
into binaries during interactions. As such, the black hole binary population in globular clusters 
will be particularly strongly affected by dynamics. 

DBHs are of great interest for gravitational wave detection since they are quite massive com- 
pared to neutron star mergers and may be detectable by the next generation of ground based de- 
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tectors at Gigaparsec distances. Strong dynamical interactions should produce many hard DBHs 
in globular clusters. It is not clear, however, that globular clusters will a strong source of DBH 
mergers since the binaries could be ejected from the cluster cores before they reach the gravita- 
tional wave regime. Sigurdsson and Hernquist (1993) |350| presented the first detailed study of 
black hole-black hole binaries in globular clusters. They used purely analytic estimates to show 
that black holes would indeed mass segregate and form a Spitzer unstable subsystem in the cluster 
core. In this subsystem black holes interact with each other to form binary and triple systems. 
Interactions between these systems and the remaining black holes remove both from the cluster 
very efficiently and Sigurdsson and Hernquist predict that there should be only zero to four black 
holes left in the cores of old galactic globular clusters, with perhaps a few more in the cluster halos. 
Remarkably, this basic picture of mass segregation, dynamical interaction, binary formation and 
ejection has not changed much since Sigurdsson and Hernquist's first publication. 

Studies of black hole binaries remained dormant until the direct N-body simulations of Portegies 
Zwart and McMillan (2000) |307) . Using smah simulations of between 2048 and 4096 particles 
with a bimodal mass distribution, 0.5-1% of bodies ten times more massive than the rest, they 
essentially confirmed the results of Sigurdsson and Hernquist. They found that the black holes 
mass segregate rapidly, form binaries through three-body interactions, and are then ejected from 
the clusters. Overall they found only ^8% of black holes were retained by their parent clusters 
while ~61% were ejected as singles and '^31% were ejected as binaries. The found, however, 
that many of the ejected, dynamically-formed black hole-black hole binaries were tight enough to 
merge in the galactic field outside the boundaries of the clusters within a Hubble time. Indeed, 
extrapolating their results to globular clusters they estimated up to one advanced LIGO detection of 
a dynamically-formed DBH binary merger per day. Thus Portegies Zwart and McMillan concluded 
that globular cluster would be an important source of DBH binaries. It must be noted, however, 
that these simulations did not include stellar evolution or natal kicks for black holes. 

During the mid-2000s the focus shifted towards the possibility that, as proposed by Miller and 
Hamilton (2002) [265', IMBHs could be formed by repeated mergers of stellar-mass black holes 
in dense cluster cores. Giiltekin et al. (2004) jl46j performed numerical three-body scattering 
experiments in a static cluster background, using the output of the last scattering experiment as the 
input of the next, to follow the merger-induced growth of black holes of various seed masses, initially 
in binaries with IOMq companions. They used the static background to calculate an encounter 
timescale between interactions and, if the quadrapole gravitational wave radiation timescale was 
less than the encounter timescale, a merger occurred. They found that it was not possible to create 
a merged black hole with a mass > 24OM0 for any reasonable seed mass. Either the proto-IMBH 
was ejected or the cluster core was totally depleted of black holes before the proto-IMBH could 
grow to a greater mass (see Figure fT3|) . They confirmed, however, that many of the ejected binaries 
would pass thorough the LISA band and then merge in the VIRGO/LIGO band within a Hubble 
time. They also found that many of the binaries would retain significant eccentricity in the LISA 
band. 

O'Leary et al. (2006) 1283) considered the same problem using a slightly different approxima- 
tion. They assumed all black holes were completely mass-segregated and formed a Spitzer- unstable 
subsystem interacting only with itself. The they followed the evolution of this subsystem, embed- 
ded in a globular cluster potential, for a range of initial conditions using a Monte Carlo selection 
of encounters and numerical integration of binary-single and binary-binary interactions. In some 
of their simulations they included a prescription for gravitational wave recoil after a DBH merger. 
Without recoil they found the could build a merged black hole of a few hundred solar masses, 
similar to Giiltekin et al. although slightly more massive. With recoil, however, the found that 
few, if any, merged products remained within the cluster and essentially no black hole growth 
occurred at all. Thus O'Leary et al. also conclude that IMBH formation through merging stellar 
mass black holes is impossible. 
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Figure 13: Total number of black holes ejected in building up to a lOOOM© IMBH as a function of 
seed mass. The four numbered curves are made assuming different cluster escape velocities. The 
dotted line is 10"^ black holes, the assumed population for a young cluster. As can be seen, a very 
large seed mass is necessary before IMBH formation by merging black holes becomes a possibility. 
Figure taken from |146j 

O'Leary et al. find, however, a large number of DBH binaries are formed through various 
interactions, many of which will merge within a Hubble time. They quantified the detection rate 
for advanced LIGO using simplified cosmological assumptions and assuming optimal orientation. 
The predict 10s of black hole merger detections per year, the exact number depending strongly on 
the initial conditions of the cluster model, and find that between 50% and 70% of these will take 
place in binaries that have been ejected from the clusters. Again, the clusters produce many black 
hole-black hole binaries but do not retain them. 

Both Giiltekin et al. and O'Leary et al. included only self-interacting black hole systems in 
their simulations. By contrast Sadowski et al. 2008 |330j modeled entire globular clusters using 
the same Monte Carlo method as Ivanova et al. |200| and assumed that the black holes remained 
in thermal equilibrium with the rest of the cluster. This is, in many ways, the opposite dynamical 
assumption to that made by O'Leary et al. since it assumes that the black holes are minimally 
mass segregated. It maximizes the number of interactions between black holes and other cluster 
stars while minimizing the number of interactions between black holes and other black holes. The 
main result of this is that there are fewer black hole-black hole binaries formed but those that do 
form are far less likely to experience a further encounter with another black hole. Interactions with 
lower-mass cluster members are less energetic and thus less likely to disrupt or eject DBHs. In 
practice the lower destruction and ejection rate wins out and Sadowski et al. predict a very large 
number of DBH binaries in globular clusters. Indeed they estimate an advanced LIGO detection 
rate of 25-3000 per year, 90% of which will occur within the boundaries of the clusters. When 
compared with the results of O'Leary et al., this proves that details of the treatment of global star 
cluster dynamics are very important. 

Two recent sets of studies have tried to address the problem of more accurate global dynamics. 
The Monte Carlo simulations of Downing et al. 2010 and 2011 [87j[88] and the direct N-body simu- 
lations of Banerjee et al. 2010 ^T\. Downing et al. use the Monte Carlo code of Giersz [12011122] to 
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44 ±6 




44 ±6 


101ow75 


18± 10 
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112 ±24 
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127 ± 16 
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175 ± 26 




184 ±29 
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64± 10 


155 ± 22 


173 ± 22 


202 ± 38 


501owl80 


103 ± 19 


205 ± 35 


294 ± 50 


453 ± 109 



Table 6: The cumulative number of BH-BH binaries after 3, 9, and 25 tj-h, and also after ITh- 
The first column gives the initial conditions. 10 or 50 in the first place indicates 10% or 50% 
primordial binaries, sol or low indicates a metallicity of z = 0.02 or z = 0.001. The number in 
the third place indicates rt/rh = 21, 37, 75 or 180. Each column is averaged over ten independent 
realizations and includes the la uncertainty. A dash in a column indicates that the cluster did not 
reach that number of half-mass relaxation times within one Hubble time. Taken from Downing et 
al. 2010 [87] 

perform 160 fully self-consistent simulations of globular clusters with 500 000 particles, a realistic 
mass function, stellar evolution, 10-50% primordial binaries, and various combinations of initial 
concentration and metallicity. The results of the simulations tend to support the approximation 
of O'Leary et al., namely the black holes mass segregate very quickly and form a dense system 
in the core of the cluster. Downing et al. find that many DBH binaries are formed, particu- 
larly by three-body interactions, but many of these binaries are either destroyed in subsequent 
interactions or are ejected from the cluster. In dense cluster DBH production is strongly time- 
dependent. Many DBHs are formed early but this number drops substantially after a few Gyrs as 
the clusters are depleted of black holes. By contrast the clusters with a low initial concentration 
do not suffer as strong a depletion in their black hole population and produce DBHs slowly but 
constantly throughout their lives. The number of black hole binaries produced is highly stochastic 
with large simulation-to-simulation variations, even between models with the same initial condi- 
tions, as shown in Table \6\ Depending on the initial conditions Downing et al. predict between 
1-100 DBH merger detections per year for advanced LIGO with an average of ~15 per year. This 
is consistent with, although slightly higher than, O'Leary et al. The detection rate is degenerate 
between some simulations with different initial conditions and may be stochastic as well. This 
means that the black hole-black hole merger detection rate alone will probably not carry much 
astrophysical information. Downing et al. also predict the possibility of eccentric DBH binaries in 
the LISA band. 

Banerjee et al. use NB0DY6 to perform direct N-body simulations of young to intermediate 
age massive clusters with between 5000 and 100 000 stars with stellar evolution and a realistic IMF. 
The also perform a series of 3000 to 4000 body simulations surrounded by a reflective boundary 
to model the core of a globular cluster. They find that the black holes mass-segregate quickly, on 
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timescales of 50-100 Myr, almost independently of the mass of the cluster. Once mass-segregated 
the black holes interact strongly, making many black hole-black hole binaries during three-body 
encounters. Like most previous authors, Banerjee et al. find that the black holes and black hole 
binaries are efficiently disrupted or ejected and that their clusters are depleted of black holes on a 
4.5 Gyr timescale. Thus they predict that current galactic globular clusters contain no black holes. 
They also estimate That intermediate age massive clusters will produce '~31 black hole binary 
merger detections per year for advanced LIGO. This is of the same order as found by Downing et 
al. and slightly higher than the rate found by O'Leary et al. while it is about an order of magnitude 
lower than that found in the previous N-body simulations of Portegies Zwart and McMillan. They 
caution, however, that this rate is not for the old globular clusters that were the focus of the 
previous studies and thus care must be taken in the comparison. 

There have been recent efforts to include Post-Newtonian (PN) dynamics into direct N-body 
simulations. This has been done both for stellar mass black holes [228 and for black holes around 
a central massive object [5]. The algorithm seems to work reasonably well but it is not fully clear 
when a PN treatment is necessary. For central massive objects where there are many stars in the 
influence radius the two-body and cluster space and time scales are not well-separated so relativistic 
dynamics may be important for the global cluster evolution. By contrast for stellar-mass objects 
the two body and cluster scales are well-separated so it is probably not important to include PN 
dynamics as part of the global simulation. 

In 2010 Ivanova et al. |198) considered the possibility of black hole-white dwarf X-ray binaries 
as UCXBs in globular clusters. This study employed scattering experiments but a rather more sim- 
plified treatment of cluster dynamics and the interaction between binaries and the cluster. Ivanova 
et al. determine that stellar collisions, exchange interactions and three-body binary formation can 
all produce black hole-white dwarf X-ray binaries, however in most cases the binary is too wide 
to become an X-ray source within a Hubble time. Thus they predict that most black hole-white 
dwarf X-ray binaries must have experienced dynamical hardening, either through fly-by encounters 
or as part of a hierarchical triple where the inner black hole- white dwarf binary has its separation 
reduced by the Kozai mechanism 222 1 . Each formation mechanism, coupled with dynamical hard- 
ening, can explain the observed number of ULXBs in globular clusters but all require that > 1% 
of black holes remain in the cluster upon formation. This places limits on the possible natal kick 
distribution of black holes. 

Unlike some of the other species mentioned, the basic behaviour of black hole-black hole binaries 
in globular clusters is well-understood. The black holes mass segregate, interact, form binaries and 
are swiftly ejected. The focus of future investigations will be on event rates and details of mass 
distributions. It seems certain, however, that globular clusters will be a major source of DBH 
binary mergers. When compared with the results for fleld DBH binaries of Belczynski et al. 2007 
[34] all of the studies cited indicate that dynamically formed binaries will dominate the DBH 
binary merger detection rate for advanced LIGO by factors of several. Some more recent results 
|30j indicate that the number of fleld black hole binaries may be larger than previously estimated 
but dynamically formed binaries are still expected to be a major source of detections. Including 
younger clusters before they have been depleted in black holes may be important for estimating 
event rates and will probably only cause them to increase. Should the next generation of ground- 
based detectors perform according to expectations, detection of a DBH binary inspiral is almost 
certain. 

5.3 Semi-empirical methods 

The small number of observed relativistic binaries can be used to infer the population of dark 
progenitor systems [38] . For example, the low- mass X-ray binary systems are bright enough that 
we see essentially all of those that are in the galactic globular cluster system. If we assume that 
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the ultracompact ones originate from detached WD-NS systems, then we can estimate the number 
of progenitor systems by looking at the time spent by the system in both phases. Let A'x be the 
number of uhracompact LMXBs and Tx be their typical lifetime. Also, let A^dot be the number of 
detached WD-NS systems that will evolve to become LMXBs, and Tdet be the time spent during 
the inspiral due to the emission of gravitational radiation until the companion white dwarf fills its 
Roche lobe. If the process is stationary, we must have 

^ = ^. (38) 

X J det 

The time spent in the inspiral phase can be found from integrating the period evolution equation 

^ = ^k,P'^/^ (39) 
at 

to obtain 

Tdet - ^ {P^^' ~ P^/') , (40) 

where 

, 96 ,8/3 G^^^ M1M2 

ko = —i2n)'—^————, (41) 
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Pq is the period at which the progenitor emerges from the common envelope, and Pc is the period 
at which RLOF from the white dwarf to the neutron star begins. Thus, the number of detached 
progenitors can be estimated from 

There are six known ultracompact LMXBs with orbital periods small enough to require a 
degenerate white dwarf companion to the neutron star. There are four other LMXBs with unknown 
orbital periods. Thus, 4 < Ax < 10. The lifetime Tx is rather uncertain, depending upon the 
nature of the mass transfer and the timing when the mass transfer would cease. A standard 
treatment of mass transfer driven by gravitational radiation alone gives an upper bound of Tx ^ 
10^ yr [316] ■ but other effects such as tidal heating or irradiation may shorten this to Tx ~ 
10^ yr [El [320]. The value of Po depends critically upon the evolution of the neutron star-main- 
sequence binary, and is very uncertain. Both fcp and Pc depend upon the masses of the white dwarf 
secondary and the neutron star primary. For a rough estimate, we take the mass of the secondary 
to be a typical He white dwarf of mass 0.4 and the mass of the primary to be 1.4 M0. Rather 
than estimate the typical period of emergence from the common envelope, we arbitrarily choose 
Pq = 2000 s. We can be certain that all progenitors have emerged from the common envelope by 
the time the orbital period is this low. The value of Pc can be determined by using the Roche lobe 
equation of Eggleton [53] : 

0.6g2/3 + In (1 + gi/3) ' 

and the radius of the white dwarf as determined by Lynden-Bell and O'Dwyer [2421 . Adopting the 
optimistic values of Ax = 10 and Tx = 10'' yr, and evaluating Equation pO]) gives Tdot ^ lO*" yr. 
Thus, we find Adct '^1- 10, which is within an order of magnitude of the numbers found through 
dynamical simulations and encounter rate estimations. 

Current production of ultracompact WD-NS binaries is more likely to arise through collisions 
of neutron stars with lower mass red giant stars near the current turn-off mass. The result of such 
a collision is a common envelope that will quickly eject the envelope of the red giant and leave 
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behind the core in an eccentric orbit. The resuh of the eccentric orbit is to hasten the inspiral of 
the degenerate core into the neutron star due to gravitational radiation |293j . Consequently, Tdot 
can be significantly shorter [201] . Adopting a value of Tdct ~ 10^ gives iVdct < 100. 

Continuing in the spirit of small number statistics, we note that there is one known radio pulsar 
in a globular cluster NS-NS binary (B2127+11C) and about 50 known radio pulsars in the globular 
cluster system as a whole (although this number may continue to grow) [238) . We may estimate 
that NS-NS binaries make up roughly 1/50 of the total number of neutron stars in the globular 
cluster system. A lower limit on the number of neutron stars comes from estimates of the total 
number of active radio pulsars in clusters, giving A^ns^ ^ 10^ [227) . Thus, we can estimate the 
total number of NS-NS binaries to be ^ 2000. Not all of these will be in compact orbits, but 
we can again estimate the number of systems in compact orbits by assuming that the systems 
gradually decay through gravitational radiation and thus 

A^compact ^compact (44) 

^NS^ ^coal cscc 

where Acompact is the number of systems in compact orbits (Porb < 2000 s), ^compact is the time 
spent as a compact system, and Tcoaiesce is the typical time for a globular cluster NS-NS binary to 
coalesce due to gravitational radiation inspiral. Adopting the coalescence time of B2127+11C as 
typical, Tcoaioscc = 2 x 10* yr [3T2j, and integrating Equation ([40|) for two 1.4 neutron stars, 
we find Acompact ^ 25. Again this value compares favorably with the values found from encounter 
rate estimations. 
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6 Prospects of Gravitational Radiation 

In the coming decades, several gravitational wave observatories will begin detecting signals from 
relativistic binaries. It is likely that several of these binaries will have their origins within globular 
clusters. In this section, we review the detectors and their expected sources from globular clusters. 

Ground-based interferometers such as the Laser Interferometer Gravitational- wave Observatory 
(LIGO) j240] in the United States, Virgo |378j in Europe, or the Large-scale Cryogenic Gravita- 
tional wave Telescope (LCGT) j230J in Japan, will be sensitive to gravitational radiation in the 
frequency range from a few Hz to a few kHz. In this frequency range, these detectors will be sensi- 
tive to the inspiral and merger of binaries containing neutron stars or black holes. These systems 
are somewhat rare in the Galactic globular cluster system, and the probability of a coalescence 
occurring within our Galaxy while the detectors are operating is vanishingly small. The enhanced 
design of LIGO was operational until the end of 2010, but did not detect any known gravitational 
wave signal. The enhanced design of Virgo will continue operating until 2012, but with the absence 
of any other interferometers operating at similar sensitivities it is unlikely that it will detect any 
known gravitational wave signals before shutting down. Both Virgo and LIGO will be upgrading 
to advanced designs which will provide significant increases in sensitivity. Like advanced LIGO 
(aLIGO) and advanced Virgo (AdV), LCGT is a next-generation detector. It is unique in that 
it will be in an underground facility with and it will operate cryogenically. All of these detectors 
should be sensitive to double neutron star mergers at distances of a few hundred Mpc, and double 
black hole mergers at distances of about a Gpc. At these distances, the number of extragalactic 
globular clusters is large enough that coalescences of neutron star or black hole binaries may be 
commonplace. Unfortunately, the angular resolution of the network of next-generation detectors 
will be much too coarse to identify any detected signal with a globular cluster. However, the grav- 
itational waveform of the coalescence of these objects depends upon the orientation of the spins of 
the components relative to the orbital angular momentum. Although it is not entirely clear if field 
binaries will tend to have their spins aligned, there is no reason for this to be the case with bina- 
ries that have been formed dynamically through exchange of compact objects. Estimates of the 
event rates of compact object binary coalescence during the advanced detector era tend to weight 
the binary neutron star coalescences toward field binaries and the binary black hole coalescences 
toward cluster binaries. Successful parameter estimation for the properties of coalescing binaries 
in advanced detectors will require accurate modeling of the waveform, including spin effects for 
both black holes and neutron stars, and tidal effects for neutron stars. 

Binary systems containing white dwarfs will be brought into contact at orbital periods of a 
few hundred seconds. At contact, the system will either coalesce on a dynamical timescale if 
the mass transfer is unstable, or it will begin to spiral out if the mass transfer is stable. In any 
case, these systems will never produce gravitational radiation in the frequency band of ground- 
based interferometers. Space-based interferometers are capable of achieving sufficient sensitivity 
in the millihertz band. Thus, the sources for space-based interferometers will include white dwarf, 
neutron star, and black hole binaries. At these orbital periods, the systems are not very relativistic 
and the gravitational waveforms can safely be described by the quadrapole formula of Peters & 
Mathews [2931 1294] . An angle averaged estimate of the typical strain amplitude is [35] 



At a typical globular cluster distance of r 10 kpc and typical chirp mass oi A4 ^ 0.5 Mq, a 
relativistic WD-WD or WD-NS binary with Porb = 400 s will have a gravitational wave amplitude 
of 10^^^. This value is within the range of proposed space-based gravitational wave observatories 
such as LISA [33] ■ 

Many globular clusters lie off the plane of the galaxy and are relatively isolated systems with 




(45) 



57 



known positions. The angular resolution of LISA improves with signal strength. By focusing 
the search for gravitational radiation using known positions of suspected sources, it is possible to 
increase the signal-to-noise ratio for the detected signal. Thus, the angular resolution of LISA for 
globular cluster sources can be on the order of the angular size of the globular cluster itself at 
/gw > 1 mHz. Consequently, the orbital period distribution of a globular cluster's population of 
relativistic binaries can be determined through observations in gravitational radiation. We will 
discuss the prospects for observing each class of relativistic binaries covered in this review. 

WD-WD binaries that are formed from a common envelope phase will be briefly visible while 
the recently revealed hot core of the secondary cools. These objects are most likely the "non- 
flickerers" of Cool et al. [62] and Edmonds et al. [92] . WD-WD binaries formed through exchange 
interactions may very well harbor white dwarfs which are too cool to be observed. In either 
case, hardening through dynamical interactions will become less likely as the orbit shrinks and 
the effective cross section of the binary becomes too small. These objects will then be effectively 
invisible in electromagnetic radiation until they are brought into contact and RLOF can begin. 
During this invisible phase, the orbital period is ground down through the emission of gravitational 
radiation until the orbital period is a few hundred seconds [351 . With a frequency of 1 to 10 mHz, 
gravitational radiation from such a binary will be in the band of LISA |39j . 

WD-NS binaries that are expected to be progenitors of the millisecond pulsars must pass 
through a phase of gravitational radiation after the degenerate core of the donor star emerges from 
the common envelope phase and before the spin-up phase begins with the onset of mass transfer 
from the white dwarf to the neutron star. The orbital period at the onset of RLOF will be on 
the order of 1 to 2 minutes and the gravitational wave signal will be received at LISA with a 
signal-to-noise of 50- 100 at a frequency of around 20 mHz for a globular cluster binary. 

Binaries with significant eccentricity will have a spectrum of harmonics of the orbital frequency, 
with the relative strength of the nth harmonic for eccentricity e given by ^293j 



where J„ is the Bessel function. The higher harmonics of sufficiently eccentric binaries (e > 0.7) 
can be detected by LISA even though the fundamental orbital frequency is well below its sensitivity 
band of 1 - 100 mHz [S^. Dynamical interactions may produce an eccentric population of 30-140 
white dwarf binaries that would be present in the LISA data after a 5 year observation [387] . 

Although the globular cluster population of NS-NS binaries is expected to be quite small 
(~ 10), they may have high eccentricities. The binary pulsar B2127-I-11C is an example of a NS- 
NS binary in a globular cluster. In terms of the unknown angle of inclination i, the companion 
mass to the pulsar is AI2 sini ~ 1 Mq and its eccentricity is e = 0.68 |238| . These binaries may also 
be detectable by LISA. If the globular cluster systems of other galaxies follow similar evolution 
as the Milky Way population, these binaries may be potential sources for LIGO as gravitational 
radiation grinds them down to coalescence. With their high eccentricities and large chirp mass, 
black hole binaries will also be good potential sources for gravitational radiation from the galactic 
globular cluster system [36l |37] . 

The relatively close proximity of the galactic globular cluster system and the separations be- 
tween individual globular clusters allows for the identification of gravitational radiation sources 
with their individual host clusters. Although the expected angular resolution of LISA is not small 
enough to allow for the identification of individual sources, knowledge of the positions of the clus- 
ters will allow for focused searches of the relativistic binary populations of the majority of the 
galactic globular clusters. Armed with a knowledge of the orbital periods of any detected bina- 
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ries, concentrated searches in electromagnetic radiation can be successful in identifying relativistic 
binaries that may have otherwise been missed. 

There may be a few double black hole systems within the Galactic globular cluster system that 
are either within the LISA band, or are sufficiently eccentric that their peak gravitational wave 
power is in a high harmonic that is within the LISA band. In either case, such systems would 
be easily detectable by a space-based gravitational wave detector due to their large chirp masses. 
However, a greater number of double black hole systems are expected to have been ejected from 
their birth globular clusters and these would now be inhabitants of the Galactic halo. If these 
systems have been formed dynamically, they are likely to be composed of two relatively massive 
black holes, and they may have chirp masses that are high enough for them to be detected at 
extragalactic distances. 
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7 Summary 



Relativistic binaries are tracers of the rich dynamical evolution of globular clusters. The properties 
of this binary population are the result of an interplay between the gravitational dynamics of large 
A''-body systems, the dynamics of mass transfer, the details of stellar evolution, and the effect of the 
gravitational field of the galaxy. The gravitational dynamics of globular clusters can enhance the 
population of short period binaries of main-sequence stars as well as inject compact objects such 
as white dwarfs and neutron stars into stellar binary systems. Once they are in such systems, the 
details of stellar evolution and mass transfer in close binary systems govern the likely end products 
of the dynamical interaction between the two stars. Furthermore, most models of the evolution of 
the core of a globular cluster rely on the gradual hardening and ejection of binary systems to delay 
the onset of core collapse. The hardening of binaries in the core of globular clusters will produce 
relativistic binaries, but it will also eventually eject these systems as they gain larger and larger 
recoil velocities in each subsequent encounter. The threshold for ejection from a globular cluster 
depends both upon the gravitational potential of the cluster itself and the gravitational potential 
of its environment generated by the Milky Way. As the globular cluster orbits the Milky Way, 
its local environment changes. Consequently, if other dynamical processes (such as gravothermal 
oscillations) do not dominate, the globular cluster's population of relativistic binaries may also 
reflect the past orbital history of the globular cluster. 

Over the last decade, observational techniques and technology have improved to the extent that 
significant discoveries are being made regularly. At this point, the bottleneck in observations of 
binary millisecond pulsars, low-mass X-ray binaries, and cataclysmic variables is time, not tech- 
nology. As these observational techniques arc brought to bear on more clusters, more discoveries 
are bound to be made. In the next decade, the possibility of using gravitational wave astronomy 
to detect relativistic binaries brings the exciting possibility of identifying the populations of elec- 
tromagnetically invisible objects such as detachcid white dwarf and neutron star binaries and black 
hole binaries in globular clusters. These observations can only help to improve the understanding 
of the complex and interesting evolution of these objects and their host globular clusters. 
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